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SUMMARY 


The  nonlinear  interaction  analysis  of  wind-driven  oblique  surface  waves  by  Lee  (J.  Fluid.  Mech., 
vol.  707,  2012,  pp.  150-190)  is  generalized  to  include  the  case  where  the  secondary-wave  frequency 
is  higher  than  the  primary  frequency.  The  nonlinear  interaction  is  initiated  by  two  three-dimensional 
fluctuations  in  wind,  which  are  induced  by  oblique  surface  waves.  If  both  wind  waves  have  the  same 
wave  speed  in  the  downwind  direction,  they  share  the  same  critical  layer.  The  primary  air  fluctuation, 
whose  steepness  of  the  order  of  the  cube  of  the  density  ratio  of  air  to  water  is  the  largest,  then  interacts 
with  a  secondary  fluctuation  of  smaller  steepness  in  the  common  critical  layer.  The  critical-layer 
interaction  generates  N  difference  modes,  where  N  is  the  smallest  integer  not  less  than  the  ratio  of 
the  secondary  to  primary  frequencies.  The  first-difference  mode  is  generated  by  a  nonlinear  coupling 
between  the  primary  and  secondary  air  fluctuations.  A  subsequent  nonlinear  interaction  between  the 
primary  fluctuation  and  the  preceding  (n—  l)th-difference  mode  generates  the  succeeding  nth-difference 
mode.  However,  the  critical-layer  interaction  between  the  primary  and  IVth-difference  modes  completes 
the  difference-mode-generation  mechanism  by  producing  a  nonlinear  effect  on  the  previous  (. N  —  l)th- 
difference  mode.  Only  the  frequencies  of  the  last  two  difference  modes,  the  (. N  —  l)th  and  IVth,  are 
lower  than  the  primary  frequency.  Their  amplitudes  can  grow  faster  than  exponentially,  and  they  soon 
approach  an  exponential-of-an-exponential  growth  stage  in  the  inviscid  and  finite-viscosity  cases. 


INTRODUCTION 

It  is  shown  by  Lee  (2012)  that  the  nonlinear  critical-layer  interactions  of  three-dimensional  fluctua¬ 
tions  in  air,  that  are  induced  by  obliquely  propagating  wind  waves,  can  significantly  enhance  the  growth 
of  one  of  the  corresponding  surface  waves.  The  critical  layer  is  a  thin  horizontal  region  surrounding  the 
critical  level  at  which  the  mean  wind  velocity  equals  the  wave  speed  in  the  downwind  direction.  The 
2  +  1  mode  interaction  in  Lee  (2012)  is  composed  of  two  wind-driven  oblique  surface  waves,  referred 
to  as  the  primary  and  secondary  waves,  and  one  difference  mode.  The  primary  wave  is  the  largest 
oblique  wave  in  the  initial  stage.  The  initial  steepness  of  a  secondary  wave  is  smaller  than  the  primary 
one.  The  2  +  1  mode  interaction  can  occur  when  both  surface  waves  have  the  same  wave  speed  in  the 
downwind  direction  and  the  frequency  of  a  secondary  wave  is  lower  than  the  primary  frequency.  The 
primary  and  secondary  waves  grow  exponentially  in  the  initial  stage  of  wind  wave  development  wherein 
the  wave  steepness  is  small  (Miles,  1957).  Once  the  steepness  of  the  primary  wave  becomes  sufficiently 
large,  the  2  +  1  mode  interaction  can  enhance  the  growth  of  the  secondary  and  difference  modes.  A 
single  difference  mode  is  generated  by  the  critical-layer  interaction  of  the  primary  and  secondary  waves 
when  the  frequency  of  the  latter  is  lower  than  the  former. 

In  this  paper,  the  case  where  the  frequency  of  a  secondary  wave  is  higher  than  the  primary  frequency 
will  be  investigated.  More  than  one  difference  mode  will  be  generated  by  critical-layer  interactions. 
Among  them,  the  difference  modes  whose  frequencies  are  lower  than  the  primary  frequency  can  be 
greatly  amplified  due  to  nonlinear  effects.  This  mechanism  may  be  responsible  for  generating  lower- 
frequency  air  fluctuations  over  surface  waves. 

A  three-dimensional  air  fluctuation  can  be  induced  by  a  surface  wave  propagating  obliquely  to 
the  downwind  direction.  The  corresponding  air  fluctuation  has  the  same  wave  characteristics,  i.e. 
frequency  and  wavenumbers,  as  the  surface  wave.  The  fluctuation  in  wind  is  governed  by  Rayleigh’s 
stability  equation  and  the  wave  motion  in  water  is  irrotational.  They  are  coupled  through  the  interface 
conditions.  As  shown  by  the  linear  theory  of  Miles  (1957),  the  wave-induced  air  fluctuation  can  be 


1 


unstable.  The  air  fluctuation  and,  consequently,  the  corresponding  surface  wave  can  grow  due  to  the 
energy  transfer  from  wind  to  air  fluctuation. 

The  linear  growth  rate  of  a  surface  wave  by  this  resonant  wind-wave  interaction  is  proportional 
to  the  density  ratio  of  air  to  water  whose  typical  value  is  O(10~3),  as  shown  by  Miles  (1957).  Since 
it  is  small,  the  leading-order  Rayleigh’s  stability  equation  becomes  singular  at  the  critical  level  where 
nonlinear  effects  first  come  into  play  (Davis  1969;  Benney  &  Bergeron  1969;  Drazin  &  Reid  2004).  If 
a  secondary  wave  has  the  same  wave  speed  in  downwind  direction  as  the  primary  wave,  both  primary 
and  secondary  fluctuations  in  air  share  the  same  critical  layer.  A  nonlinear  interaction  between  them 
can  occur  when  the  steepness  of  a  primary  wind  wave  becomes  of  the  order  of  the  cube  of  the  density 
ratio  of  air  to  water  as  found  in  Lee  (2012).  The  secondary  wave  can  be  any  oblique  wind  wave  of  the 
same  wave  speed  in  downstream  direction  whose  steepness  is  smaller  than  the  primary  one.  By  using 
a  matched  asymptotic  analysis  as  in  Reutov  (1980),  Goldstein  Sz  Choi  (1989),  Goldstein  &  Lee  (1992), 
Wu  (1992),  Wu  &  Stewart  (1996),  Lee  &;  Wundrow  (2011),  and  others,  Lee  (2012)  showed  that  the 
2  +  1  mode  critical-layer  interaction  induced  by  two  oblique  wind  waves  could  enhance  the  growth  of 
a  secondary  wind  wave. 

In  a  2  +  1  mode  interaction,  the  primary-wave  frequency  is  higher  than  the  secondary  frequency. 
The  critical-layer  interaction  between  the  primary  and  secondary  fluctuations  generates  a  difference 
mode  whose  frequency  and  wavenumbers  are  equal  to  the  differences  between  the  primary  and  sec¬ 
ondary  values.  Its  amplitude  becomes  as  large  as  the  secondary  amplitude  in  air,  but  smaller  than 
the  latter  in  water  since  the  difference  mode  does  not  satisfy  the  free-surface-wave  dispersion  relation. 
In  addition,  the  nonlinear  coupling  between  the  primary  and  difference  modes  produces  a  parametric- 
growth  effect  on  the  secondary  wind  wave.  The  primary  wave  remains  linear  throughout  the  whole 
stage.  However,  the  growth  rate  of  a  secondary  wave  becomes  much  larger  than  its  linear  growth  rate, 
and  both  secondary  and  difference  mode  amplitudes  can  grow  like  the  exponential-of-an-exponential. 
The  dynamics  of  a  2+1  mode  interaction  are  similar  to  those  of  the  phase-locked  interaction  of  Wu  & 
Stewart  (1996). 

In  the  present  study,  the  critical-layer  analysis  of  Lee  (2012)  is  extended  so  that  the  secondary 
frequency  can  be  higher  than  the  primary  frequency.  The  nonlinear  interaction  is  still  induced  by  two 
wind-driven  oblique  surface  waves.  As  in  Lee  (2012),  the  magnitude  of  a  primary  oblique  wave  is  of 
the  order  of  the  cube  of  the  density  ratio  of  air  to  water,  and  that  of  a  secondary  oblique  wave  is 
smaller  than  the  primary  one.  The  critical-layer  interaction  generates  N  difference  modes,  where  N  is 
the  smallest  integer  not  less  than  the  ratio  of  the  secondary  to  primary  frequencies.  This  interaction  is 
composed  of  two  wind-driven  surface  waves  and  N  nonlinearly  generated  difference  modes  and  might 
be  termed  the  ‘2  +  N  mode  critical-layer  interaction.’  The  2  +  1  mode  interaction  in  Lee  (2012)  is  the 
simplest  case  when  N  =  1  that  involves  only  three  modes.  The  analysis  of  a  generalized  2  +  N  mode 
interaction  shows  that  both  primary  and  secondary  waves  remain  linear  if  N  >  2.  The  first  N  —  2 
difference  modes  also  grow  exponentially  although  they  are  nonlinearly  generated.  The  amplitudes  of 
the  last  two  difference  modes,  the  (N  —  l)th  and  Nth,  are  determined  by  coupled  integral  equations. 
They  can  grow  faster  than  exponentially.  Numerical  results  will  show  that  they  can  grow  like  the 
exponential-of-an-exponential,  similar  to  the  secondary  and  difference  modes  in  the  N  =  1  case.  The 
exponential-of-an-exponential  growth  only  occurs  to  the  modes  whose  frequencies  are  lower  than  the 
primary-wave  frequency.  They  are  the  (N  —  l)th  and  Nth  difference  modes  if  N  >  2,  or  the  secondary 
and  difference  modes  when  N  =  1. 

The  linear  and  nonlinear  effects  appear  at  the  same  order  in  the  amplitude  equations  in  these 
analyses.  Thus,  the  analysis  covers  both  linear  and  nonlinear  stages  and  it  describes  how  the  nonlinear 
critical-layer  interaction  gradually  evolves  from  the  interaction  between  linear  wind  waves  of  small 
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steepnesses.  In  the  initial  stage  where  the  primary- wave  steepness  is  smaller  than  0(the  density  ratio 
cubed),  the  nonlinear  amplitude  equations  reduce  to  the  linear  equations. 

The  2  +  IV  mode  analysis  (N  >  1)  shows  that  a  nonlinear  interaction  can  occur  even  when  the 
steepness  of  the  largest  primary  oblique  wave  is  relatively  small,  i.e.  0(the  density  ratio  cubed) 
or  O(10~9).  The  nonlinear  critical-layer  interaction  could  be  an  efficient  mechanism  by  which  the 
difference  modes  and  secondary  wave  (when  N  =  1)  quickly  grow  as  large  as  the  primary  wave. 
Nonlinear  interactions  in  subsequent  stages  will  be  more  involved  since  all  modes  are  of  the  same 
order  of  magnitude.  Lee  &  Wundrow  (2011)  show  that  the  modal  growth  can  be  governed  by  a  cubic 
self- interaction  mechanism  (Goldstein  &  Choi  1989)  if  the  wave  steepness  becomes  0(the  density  ratio 
raised  to  the  power  5/2). 

The  findings  of  enhanced  growth  of  oblique  waves  are  consistent  with  previous  instability  wave 
analyses.  It  has  been  shown  by  various  studies  that  nonlinear  critical-layer  interactions  of  oblique 
instability  waves  in  single-fluid  flows  produce  much  larger  growth  effect  than  the  linear  interaction 
does.  Meanwhile,  the  critical-layer  analysis  of  a  plane  surface  wave,  that  propagates  along  the  downwind 
direction,  by  Reutov  (1980)  shows  that  the  nonlinear  effect  reduces  the  growth  rate,  even  when  the 
wave  steepness  is  as  small  as  0(the  density  ratio  squared)  (also  in  Alexakis,  Young  &  Rosner  2004). 

It  is  also  found  in  Lee  (2012)  that  a  linear  plane  wave  remains  passive  although  it  is  larger  than 
the  primary  oblique  wind  wave  of  0(the  density  ratio  cubed)  in  the  2  +  1  mode  interaction  stage.  The 
direct  critical-layer  coupling  between  the  plane  and  oblique  wind  waves  does  not  produce  a  nonlinear 
effect  outside  the  critical  layer  since  the  frequency  of  the  plane  wave  is  always  lower  than  the  frequencies 
of  oblique  waves  of  an  equal  wave  speed  in  downwind  direction.  However,  the  plane  wind  wave  can 
interact  with  the  nonlinearly  generated  difference  modes  whose  frequencies  are  lower  than  the  plane- 
wave  frequency  through  the  phase-locked  interaction  of  Wu  Sz  Stewart  (1996). 

The  analysis  in  this  paper  is  similar  to  that  in  Lee  (2012).  In  fact,  the  previous  2  +  1  mode  analysis 
can  be  included  in  the  current  2  +  N  mode  analysis  as  a  special  case  where  N  =  1.  For  simplicity, 
however,  the  presentation  in  the  main  body  will  be  given  for  N  >  2.  Parts  of  the  analysis  and  discussion 
that  were  presented  in  Lee  (2012)  will  not  be  repeated.  We  have  tried  to  keep  the  same  notations  as 
in  Lee  (2012)  for  easy  cross-referencing.  However,  a  change  occurred  for  the  normalization  of  variables 
in  the  final  amplitude  equations.  In  Lee  (2012),  the  variables  in  the  final  results  are  normalized  by 
using  the  properties  of  the  primary  and  secondary  waves  of  each  set  of  interactions,  which  simplified 
the  amplitude  equations.  In  this  paper,  the  properties  of  a  wind  wave  with  a  30°  propagation  angle  are 
used  in  all  cases  so  that  the  numerical  results  of  different  sets  of  interactions  can  be  directly  compared. 
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FORMULATION 


A  right-handed  Cartesian  coordinate  system  (x,y,z)  with  unit  vectors  (i,j,k)  is  attached  to  the 
mean  water  surface  with  x  in  the  downwind  or  streamwise  direction  of  the  wind,  y  in  the  spanwise 
lateral  direction,  and  z  in  the  vertical  direction,  as  shown  in  Figure  1.  We  suppose  that  the  mean  wind 
is  two-dimensional  and  there  is  no  mean  motion  in  the  water.  All  lengths,  velocities,  pressures  and 
time  in  both  fluids  are  non-dimensionalized  by  Aa,  Ua,  paU %  and  A a/Ua,  respectively,  where  Aa  and 
Ua  will  be  chosen  as  in  equation  (8).  The  Reynolds  numbers  are  assumed  to  be  large,  so  viscous  effects 
only  enter  through  the  critical  layer.  The  effect  of  surface  tension  is  not  included.  The  air  and  water 
are  assumed  incompressible  and  their  density  ratio  is  defined  as 

Pa/Pw  =  cr ,  (1) 

where  the  subscripts  a  and  w  denote  the  quantities  of  air  and  water,  respectively,  and  a  <C  1  charac¬ 
terizes  the  small  density  ratio  of  air  to  water. 

Perturbation  to  the  air-water  interface  about  the  mean  water  surface,  z  =  0,  is  described  by 

z  =  eh(x,y,t),  (2) 

where  e  is  a  measure  of  a  wave  steepness  and  h  is  a  scaled  wave  elevation.  The  scaling  relation  between 
e  and  a  will  be  determined  by  the  dynamics  in  the  critical  layer.  The  critical  layer  surrounds  the 
critical  level  where  the  mean  wind  velocity  equals  the  wave  speed  in  the  streamwise  direction. 

The  total  velocity  u  =  iu  +  jv  +  kw  and  pressure  p  in  air  are  composed  of  a  steady  wind  and  mean 
pressure  plus  wave-induced  perturbations  u  =  iu  +  jv  +  kw  and  p, 

u  =  iU(z)  +  eu(x,t),  p  =  Pq  —  Gz  +  ep  (x,t),  (3) 

where  U (z)  is  the  steady  wind  and  Pq  denotes  the  mean  pressure.  The  gravity  parameter  is  defined  as 

G^gAa/Ul  (4) 

where  g  denotes  the  acceleration  due  to  gravity.  It  is  assumed  that  the  mean  velocity  at  the  air-water 
interface  vanishes,  i.e. 

u(  0)  =  o.  (5) 

The  fluid  motion  in  water  is  irrotational.  The  velocity  and  pressure  where  z  <  eh  are  expressed  as 

u  =  eV<j>(x,t)i  p  =  P0  -  {l/a)Gz  +  {e/a)pw.  (6) 

Motions  of  air  and  water  are  coupled  at  the  interface  defined  by  (2).  The  kinematic  interface  condition 
requires  that  the  interface  be  a  material  surface  for  both  fluid  flows.  The  pressure  at  the  interface  must 
be  continuous  according  to  the  dynamic  interface  condition. 

A  logarithmic  profile  is  used  for  the  mean  wind  as  in  Miles  (1957,  1993),  Morland  &  Saffman  (1993), 
Janssen  (2004),  Lee  (2012),  and  others, 


U  =  2.51n(9.025z  +  1).  (7) 

The  velocity  and  length  in  the  above  turbulent  wind  profile  are  non-dimensionalized  by 

Ua  =  u *,  Aa  =  9.025zo,  (8) 
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Figure  1.  (a)  Schematic  diagram  of  a  flow  field,  (b)  Wave  field  of  a  2  +  2  mode  interaction  with  wave 
propagation  angles  of  9P  =  30°,  9S  =  —57°,  9d\  ~  —79°  and  9^2  ~  83.5°,  for  the  primary,  secondary, 
lst-difference  and  2nd-difference  modes.  The  wave  crests  are  plotted  as  solid  lines  and  troughs  are  as 
dashed  lines. 


where  u*  denotes  the  friction  velocity  and  za  is  the  roughness  length.  It  is  worth  mentioning  that  the 
present  analysis  is  developed  for  general  wind  profiles.  The  specific  form  (7)  will  only  be  used  when 
the  final  equations  are  numerically  evaluated. 
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NONLINEAR  CRITICAL-LAYER  INTERACTIONS 


A  nonlinear  interaction  mechanism  that  is  termed  the  ‘2  +  N  mode  critical-layer  interaction’  will 
be  described  in  this  section  before  we  show  detailed  analyses.  A  2  +  N  mode  interaction  is  composed 
of  two  wind-driven  oblique  surface  waves  and  N  nonlinear ly  generated  difference  modes.  Diagrams  of 
interactions  for  various  values  of  N  are  shown  in  Figures  2-4.  A  representative  wave  field  by  a  2  +  2 
mode  interaction  is  illustrated  in  Figure  1(b). 

For  each  wind-driven  surface  wave,  there  exists  a  corresponding  air  fluctuation  of  the  same  frequency 
and  wavenumbers.  The  air  fluctuation  that  is  not  associated  with  a  wind  wave  can  also  produce 
a  perturbation  in  water.  The  motion  of  air  fluctuation,  that  is  governed  by  Rayleigh’s  equation,  is 
coupled  with  the  irrotational  motion  in  water  through  the  interface  conditions  on  the  water  surface 
(Miles  1957).  During  the  nonlinear  stage  of  this  analysis,  only  the  dynamics  of  air  fluctuations  inside 
the  thin  critical  layer  are  governed  by  nonlinear  equations.  Outside  the  critical  layer,  the  nonlinear 
terms  in  the  momentum  equations  become  insignificant  and  the  fluid  motions  are  described  by  linear 
equations. 

The  dispersion  relations  of  wind  waves  are  affected  by  the  dynamics  of  fluctuations  in  the  air.  This 
effect  becomes  of  a  higher  order  in  the  present  small  density  ratio  case.  The  leading-order  dispersion 
relations  of  wind  waves  equal  the  dispersion  relations  of  free-surface  waves.  However,  the  wavenumbers 
and  frequencies  of  nonlinearly  generated  modes  are  determined  by  nonlinear  interactions.  If  they  do  not 
satisfy  the  free-surface- wave  dispersion  relations  (at  the  leading  order),  the  corresponding  perturbations 
in  water  will  disappear  when  the  nonlinearly  generated  air  fluctuations  or  nonlinear  critical-layer  effects 
vanish. 

It  is  assumed  that  a  primary  wind  wave,  denoted  by  the  letter  P  in  Figures  2-4,  is  the  largest 
oblique  surface  wave  in  an  initial  linear  wave  field.  It  is  periodic  in  the  form 

exp  [icq,  (a;  -  ct)  +  i/3py\,  (9) 

where  i  =  \J—  1,  c  is  the  wave  speed  in  the  streamwise  direction,  and  ap  and  f3p  are  the  streamwise  and 
spanwise  wavenumbers,  respectively.  A  secondary  wind  wave  of  the  same  c,  denoted  by  S  or  S  in  the 
figures,  is  proportional  to 

exp[ias(x  -  ct)  +  i/3sy\.  (10) 

The  streamwise  wavenumber  as ,  or  frequency  ujs(=  asc),  of  a  secondary  wave  can  be  arbitrarily  chosen. 

Both  primary  and  secondary  waves  satisfy  the  free-surface-wave  dispersion  relation  at  the  leading 
order.  The  relation  given  in  Lee  (2012)  is  written  as 

/32  =  [c4a2/{G2  tanh2(7?Ff)}  —  l]  a2 ,  or  c2  =  (G^/a? )  tanh^Lf),  (11) 

for  i  =  {p,  s},  where 

h  =  +  P?,  (12) 

H  is  the  depth  of  water,  and  G  is  defined  by  (4).  The  quantities  of  the  primary  and  secondary  waves 
are  denoted  with  the  subscripts  p  and  s,  respectively.  The  wave  propagation  angle  with  respect  to  the 
streamwise  direction,  9t ,  becomes 

9t  =  arctan(/3j/ at),  cos  9t  =  ar/ sin <9*  =  /Vt*,  for  *  =  { P ,  «},  (13) 

where  —ir/2  <  9t  <  7r/2. 
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(14) 


A  contour  line  of  constant  c  from  (11)  is  plotted  in  Figures  2(a,  b)-4(a,  b)  when 

c  =  10,  G  =  0.05,  H  =  1035. 

For  all  wind  waves  on  this  curve,  c  =  10.  Note  that  d\j3z\/dat  >  0  and  d\9l\/dal  >  0  for  i  =  {p,  s}, 
which  can  be  proved  from  (11)-  (13).  All  numerical  results  are  evaluated  using  the  parameter  values 
listed  in  (14).  The  large  value  of  H  produces  the  results  in  deep  water.  The  contour  line  is  replotted  in 
Figures  2(c-h)-4(c-h)  as  a  dotted  curve.  The  primary-wave  propagation  angles,  9P,  are  30°,  50°  and 
70°  in  Figures  2,  3  and  4,  respectively. 

Because  the  secondary  wave  has  the  same  c  as  the  primary  wave,  they  share  the  same  critical  level. 
As  found  in  Lee  (2012),  a  nonlinear  interaction  between  these  initially  linear  wind  waves  first  occurs 
in  their  common  critical  layer  when  the  primary-wave  steepness,  e,  becomes  0(o"3).  The  steepness  of 
the  secondary  wave  is  assumed  to  be  smaller  than  the  primary  one. 

All  harmonics,  sum  and  difference  modes  of  the  primary  and  secondary  waves,  or  of  any  two 
interacting  modes,  are  generated  by  quadratic  interactions  inside  the  critical  layer  (details  will  be 
given  later  in  (76)).  Among  them,  only  the  difference  mode,  whose  frequency  and  wavenumbers  are 
equal  to  the  differences  between  those  of  two  interacting  modes,  becomes  finite  at  the  outer  edges  of 
the  critical  layer  (as  will  be  shown  in  (83)).  The  other  nonlinearly  generated  modes  become  negligibly 
small  outside  the  critical  layer.  The  difference-mode  amplitude  becomes  as  large  as  the  secondary- 
wave  amplitude  in  the  air.  If  the  difference  mode  is  an  unsteady  oblique  mode,  it  can  interact  with  the 
primary  wave.  As  will  be  described  in  the  following  subsections,  the  number  of  difference  modes,  N, 
generated  by  a  2  +  N  mode  interaction  becomes 

N  =  ceiling  (us/up)  =  ceiling^ /ap),  (15) 

where  the  ceiling  function  denotes  the  smallest  integer  not  less  than  its  argument,  and  ujp  =  apc.  In 
other  words,  the  critical-layer  interaction  generates  N  difference  modes,  if 

(N—1)ojp<lus<Nlop,  (N—l)ap<as<Nap,  (16) 

where  N  >  1.  All  nonlinearly  generated  difference  modes  have  the  same  c  as  the  primary  and  secondary 
waves.  Since  a  plane  wind  wave  does  not  directly  interact  with  oblique  wind  waves  in  the  nonlinear 
stage  of  this  analysis,  ap  and  as  must  be  greater  than  the  wavenumber  of  a  corresponding  plane  wave, 

«2 d  =  (G/c2)  tanh(c>!2D^f)- 

2  +  1  mode  critical-layer  interaction:  enhanced  growth  of  a  secondary  wave 

The  2  +  1  mode  interaction  between  two  wind  waves  and  a  single  nonlinear  mode  can  occur  when 

N  =  1  or  as  <  ap,  (17) 

as  studied  in  Lee  (2012).  The  nonlinear  interaction  between  the  primary  and  secondary  modes,  denoted 
by  P  and  S  in  Figures  2(c,  d)-4(c,  d),  generates  a  difference  mode,  denoted  by  D.  Its  wavenumbers 
are  ( ap  —  as,/3p  —  (3S )  and  velocity  fluctuation  is  periodic  in  the  form 

exp[i(ap  -  as)(x  -  ct )  +  i(/3p  -  f3s)y\,  (18) 

and  the  propagation  angle  becomes 

9d  =  arctan[(/3p  -  /3s)/(ap  -  as)].  (19) 
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Figure  2.  Diagrams  of  2  +  IV  mode  critical-layer  interactions  when  9p  =  30°.  (a,  b)  The  contour  line 
of  the  free-surface-wave  dispersion  relation  of  constant  c  =  10  in  the  ( a ,  (3)  and  (a,  9)  coordinates, 
with  G  =  0.05  and  H  =  1035.  Diagrams  of  (c,  d)  2  +  1;  (e,  f)  2  +  2;  (g,  h)  2  +  3  mode  critical-layer 
interactions.  The  a  represents  the  streamwise  wavenumbers,  ap,  as  and  |as  —  nap |,  for  the  primary, 
secondary  and  nth-difference  modes,  and  f3  and  8  are  the  spanwise  wavenumber  and  wave  propagation 
angle  of  the  corresponding  mode. 
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iagrams  of  2  +  N  mode  critical- layer  interactions  when  6p  =  50°.  (a,  b)  The  same  c  =  10 
as  in  Figures  2(a,  b).  Diagrams  of  (c,  d)  2  +  1;  (e,  f)  2  +  2;  (g,  h)  2  +  3  mode  critical-layer 
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Figure  4.  Diagrams  of  2  +  N  mode  critical- layer  interactions  when  9p  =  70°.  (a,  b)  The  same  c  =  10 
contour  line  as  in  Figures  2(a,  b).  Diagrams  of  (c,  d)  2  +  1;  (e,  f)  2  +  2;  (g,  h)  2  +  3  mode  critical-layer 
interactions. 
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Figure  5.  Conditions  on  (9p,  9S )  for  various  2  +  N  mode  interactions.  In  regions  (iQi&),  N  =  1;  (iia,b), 
IV  =  2;  (iiia,b),  iV  =  3;  (iva]ft),  IV  =  4.  The  regions  where  N  =  5,  6  and  7  are  plotted  without  labels. 
Border  lines  for  N  >  8  are  not  plotted  since  regions  become  too  narrow.  Parameter  values  in  (14)  were 
used. 


Meanwhile,  the  critical-layer  interaction  between  the  primary  and  difference  modes  produces  a  non¬ 
linear  effect  on  the  secondary  wind  wave.  The  amplitudes  of  the  secondary  and  difference  modes  are 
governed  by  an  integro-differential  and  integral  equations,  respectively,  and  they  can  grow  faster  than 
exponentially.  The  primary  wave  remains  linear.  The  frequencies  of  the  secondary  and  difference 
modes,  us  and  a id  (=  up  —  ws),  are  lower  than  the  primary  frequency  up. 

The  primary  wave  P  in  Figures  2(c,  d)-4(c,  d)  can  also  interact  with  any  secondary  wave  on  blue 
curves  to  generate  a  difference  mode  on  red  curves.  Any  combination  of  (9P,  9S)  in  both  triangular 
regions  (ia  &)  in  Figure  5  produces  a  2  +  1  mode  interaction.  The  contour  lines  of  constant  9 j  are  shown 
in  regions  (ia,b)  in  Figure  6(a).  For  simplicity,  we  only  show  the  results  when  9P  >  0  in  the  figures. 

2  +  2  mode  critical-layer  interaction:  generation  of  two  difference  modes 

More  than  one  difference  mode  is  generated  if  the  secondary  frequency  is  larger  than  the  primary 
frequency.  If 

N  =  2  or  ap  <  as  <  2 ap,  (20) 

the  critical-layer  interaction  between  the  primary  and  secondary  modes,  P  and  S  in  Figures  2(e,  f)-4(e, 
f),  generates  a  first-difference  mode  D\ .  Its  wavenumbers  are  (as  —  ap,/3s  —  j3p)  and  it  is  periodic  in 
the  form  exp[i(as  —  ap)(x  —  ct)  +  i(j3s  —  f3p)y\,  and  the  propagation  angle  9^  equals  9 ^  given  by  (19). 
Then,  the  subsequent  quadratic  coupling  between  P  and  D\  produces  a  second-difference  mode  +>2. 
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Figure  6.  Contour  lines  of  constant  wave  propagation  angles  of  difference  modes  in  the  (Qp,  0S)  plane, 
(a)  9di ;  (b)  9  42'-  (c)  9d 3-  Border  lines  of  regions  (ia,fe),  (iia,&)  and  (iiia,b)  shown  in  Figure  5  are  replotted 
as  black  dotted  curves.  Regions  (iiia,&)  are  located  above  and  below  regions  (iia,b)  as  in  Figure  5.  Units 
of  contour-line  labels  are  degrees. 
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With  wavenumbers  of  (2 ap  —  aSl  2(3P  —  (3S ),  it  is  proportional  to 

exp [i( 2 -  as)(x  -  ct )  +  i(2 /3P  -  / 3s)y }.  (21) 

The  propagation  angle  of  D-2  becomes 

Od 2  =  arctan[(2/3p  -  f}s)/ (2 ap  -  as)].  (22) 

Since  the  wavenumber  differences  between  P  and  D2  are  equal  to  the  wavenumbers  of  D\,  the  next 
nonlinear  interaction  of  the  primary  and  2nd-difference  modes  produces  a  nonlinear  effect  on  the  lst- 
difference  mode.  Both  primary  and  secondary  waves  remain  linear.  The  amplitudes  of  the  first  and 
second  difference  modes  are  governed  by  integral  equations,  and  they  can  grow  faster  than  exponentially 
as  will  be  shown  later.  The  frequencies  of  both  difference  modes,  u>di  (=  cos  —  top)  and  Ud2  (=  2 ujp  —  ujs), 
are  lower  than  the  primary- wave  frequency  uip.  The  wave  field  by  the  2  +  2  mode  interaction  of  four 
modes  P ,  S,  D\  and  D2  in  Figures  2(e,  f)  is  illustrated  in  Figure  1(b). 

The  2  +  2  mode  interaction  by  P  and  S  in  Figures  2(e,  f)  generates  the  difference  modes  D\  and 
D2.  The  interaction  of  P  and  any  secondary  wave  located  on  green  solid  or  dashed  curve  in  Figures 
2(e,  f)-4(e,  f)  generates  the  1st  and  2nd  difference  modes  on  the  corresponding  blue  and  red  curves. 
The  2  +  2  mode  interaction  occurs  in  regions  (iia,&)  in  Figure  5.  The  contour  lines  of  constant  6d\  and 
0d2  are  plotted  in  Figures  6(a,  b). 

It  is  interesting  to  note  that  0d2  can  become  zero  in  region  (iia)  in  Figure  6(b).  Figures  2(f)— 4(f) 
also  show  that  a  2  +  2  mode  interaction  can  generate  a  plane  second-difference  mode.  The  condition 
for  this  case  will  be  given  later. 

The  2nd-difference  mode  can  be  a  wind  wave  as  shown  in  Figures  3(e,  f)  and  4(e,  f)  for  the  9P  =  50° 
and  70°  cases.  The  2  +  2  mode  interaction  by  P  and  S  generates  the  difference  modes  D\  and  D2. 
Since  the  2nd-difference  mode  D2  is  located  on  the  black  dotted  curve,  it  satisfies  the  free-surface-wave 
dispersion  relation  at  the  leading  order.  This  is  an  interesting  case  where  a  3  +  1  mode  interaction 
between  three  wind  waves  and  one  difference  mode  occurs. 

2  +  N  mode  critical-layer  interaction:  generation  of  N  difference  modes 

Previous  2  +  2  mode  interaction  mechanism  can  be  extended  for  an  arbitrary  value  of  N  >  2.  Two 
wind  waves  and  N  difference  modes  are  involved  in  a  2  +  N  mode  interaction.  Diagrams  of  interactions 
when  N  =  3  are  shown  in  Figures  2(g,  h)-4(g,  h). 

When 

N  >2  and  (N—  l)ctp  <  as  <  Nap,  (23) 

the  lst-difference  mode  D\  in  the  Figures  is  generated  by  the  nonlinear  interaction  of  the  primary  and 
secondary  modes,  P  and  S',  as  in  the  2  +  2  mode  interaction.  The  subsequent  quadratic  interaction 
between  the  primary  wave  P  and  the  preceding  ( n  —  l)th-difference  mode  Dn_  1  keeps  generating  the 
succeeding  ?rth-difference  mode  Dn  until  the  last  IVth-difference  mode  Djy  is  generated.  For  1  <  n  < 
N  —  1,  the  wavenumbers  of  the  nth-difference  mode  become  (as  —  nap, /3S  —  n/3p)  and  the  velocity 
fluctuation  is  periodic  in  the  form 

exp[i(as  —  nap){x  —  ct)  +  i ((3S  —  n/3p)y\  if  1  <  n  <  N  —  1,  (24) 

Its  wave  propagation  angle  becomes 

@dn  =  arctan[((Ss  —  nj3p)/{as  —  nap)\  if  1  <  n  <  IV  —  1.  (25) 
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(26) 


The  ATh-difference  mode,  with  wavenumbers  of  ( Nap  —  as,  N/3P  —  f3s ),  is  periodic  in  the  form 

exp[i(iVap  —  as)(x  —  ct)  +  i  (N  /3P  —  (3s)y\  if  n  =  N, 

and  its  wave  propagation  angle  becomes 

0dN  =  arctan[(AI/lp  -  /3s)/(Nap  -  as)].  (27) 

As  expected  from  the  2  +  1  and  2  +  2  mode  interactions,  the  difference-mode-generation  mechanism 
is  completed  by  the  interaction  of  the  primary  and  A+h-difference  modes  since  their  wavenumber  dif¬ 
ferences  are  equal  to  the  wavenumbers  of  the  (N  —  l)th-difference  mode.  Thus,  the  coupling  between 
P  and  Djy  produces  a  nonlinear  effect  on  Djv-i-  The  linear  primary  and  secondary  waves  grow  ex¬ 
ponentially.  Although  it  is  nonlinearly  generated,  the  nth-difference  mode  also  grows  exponentially 
for  1  <  n  <  N  —  2.  Their  frequencies,  ujdn  (=  <^d(n- 1)  —  wP),  are  higher  than  the  primary  frequency. 
However,  the  frequencies  of  the  (AT  —  l)th  and  A+h  difference  modes,  c ^d(N-i)  (=  ws  —  (N  —  l)up)  and  u )dN 
(=  Noop  —ujs),  are  lower  than  oop.  It  will  be  shown  later  that  the  amplitudes  of  these  last  two  difference 
modes  are  governed  by  coupled  integral  equations  and  they  can  grow  faster  than  exponentially. 

In  regions  (iiia,fe)  in  Figure  5  where  2 ap  <  as  <  3 ap,  three  difference  modes  are  generated  by  a  2  +  3 
mode  interaction.  A  2  +  4  mode  interaction  can  occur  in  regions  (i va,b)  where  3ap  <  as  <  4ap.  For 
2  +  3  mode  interactions,  the  contour  lines  of  constant  0d\ ,  &d2  and  #cZ3  are  plotted  in  regions  (hia,b)  in 
Figure  6.  Regions  (iiia,fe)  are  located  above  and  below  regions  (iia,b)  as  in  Figure  5. 

The  present  analysis  will  show  that  the  faster-than-exponential  growth  only  occurs  to  two  modes 
whose  frequencies  are  lower  than  the  primary  frequency.  They  are  the  secondary  and  difference  modes 
in  a  2  +  1  mode  interaction,  and  the  (A"  —  l)th  and  Nth  difference  modes  in  2  +  N  mode  interactions 
with  N  >2.  When  N  >  2,  the  first  N  —  1  quadratic  interactions  between  the  primary  and  secondary/d¬ 
ifference  modes  generates  A"  —  1  difference  modes  whose  frequencies  are  successively  decreased  by  the 
amount  of  the  primary  frequency.  The  frequency  of  the  (N  —  l)th-difference  mode  first  becomes  lower 
than  the  primary  frequency.  Since  the  frequency  of  the  ATh-difference  mode  equals  the  primary  fre¬ 
quency  minus  the  ( N  —  l)th-difference  frequency,  it  is  also  lower  than  the  primary  frequency.  The 
(N  —  l)th  and  Nth  difference  modes  are  analogous  to  the  secondary  and  difference  modes,  respec¬ 
tively,  of  a  2  +  1  mode  interaction.  The  critical-layer  interaction  between  the  primary  wind  wave  and 
these  two  modes  enables  the  latter  to  grow  faster  than  exponentially,  and  eventually  to  grow  like  the 
exponential-of-an-exponential.  The  critical-layer  analysis  and  resulting  amplitude  equations  for  N  >2 
will  be  presented  in  the  following  sections. 

As  shown  in  Appendix  F.2  in  Lee  (2012),  the  direct  interaction  between  the  plane  and  oblique 
wind  waves  does  not  produce  nonlinear  effects  outside  the  critical  layer.  However,  the  phase-locked 
interaction  (Wu  &  Stewart  1996)  can  occur  between  the  plane  wave  and  nonlinearly  generated  difference 
modes  whose  frequencies  are  lower  than  the  plane-wave  frequency.  Additional  difference  modes  are 
generated  by  critical-layer  interactions  between  the  plane  and  original  difference  modes.  The  frequencies 
of  the  original  and  additional  difference  modes  are  lower  than  the  plane  wave  frequency. 
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UNSTEADY  FLOW  OUTSIDE  THE  CRITICAL  LAYER 


The  unsteady  flow  in  the  outer  region  is  composed  of  a  primary  and  secondary  wind  wave  and  N 
difference  modes  where  N  is  determined  by  (15).  The  difference  modes  are  generated  by  nonlinear 
critical-layer  interactions.  All  difference  modes  as  well  as  the  primary  and  secondary  waves  have  the 
same  wave  speed  in  the  streamwise  direction,  c.  The  difference  modes  are  assumed  not  to  satisfy  the 
free-surface-wave  dispersion  relation  (11).  There  are  special  cases  where  one  of  them  satisfies  (11)  as 
mentioned  before. 

The  analysis  in  this  section  closely  follows  that  in  §3  in  Lee  (2012).  One  important  difference  is 
that  the  secondary  wave  is  linear  when  N  >  2,  although  it  was  nonlinear  in  Lee  (2012)  when  N  =  1. 
The  present  analysis  also  includes  multiple  difference  modes.  The  primary  wave  is  linear  in  both  cases. 


Unsteady  flow  in  the  water 

The  velocity  potential  <f>  and  pressure  perturbation  pw  in  the  water,  and  the  interface  surface  elevation 
h  can  be  written  as, 


(  4>  \ 

Pw 

\'h  ) 


=  Br 


(  ,  I®  \ 
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V  J 
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£s  +  a  5  ^  Bdr 
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^  <t>dn  ^ 
Pdn 

V  J 


£dn  +  C.C.  +  .  •  •  ,  (28) 


where  c.c.  denotes  the  complex  conjugate,  and 


£%  =  exp[i(aiX+/3ty)\  for  i  =  {p,s,dN}, 

£dn  =  exp[i(-adnA-/3dny)]  if  l<n<AT-l, 

adn  =  nap  -  as,  /3dn  =  n/3p  -  f3s,  for  1  <  n, 
X  =  x  —  ct, 


(29) 

(30) 

(31) 


with  X  denoting  the  normalized  streamwise  coordinate  in  a  reference  frame  moving  with  the  common 
c. 

The  subscripts  p.  s  and  dn  denote  the  primary,  secondary  and  ?rth-difference  modes,  respectively. 
The  streamwise  wavenumbers,  ap,  as,  —adl,  ...,  —ad(N-i),  and  adN  are  positive.  But,  the  spanwise 
wavenumbers  fdp ,  f3s  and  j3dn  can  be  positive  or  negative  as  illustrated  in  Figures  2-4.  As  in  (12)  and 

(13), 

7 dn  =  \J ajn  +  P2dn  for  1  <  n,  (32) 

@dn  arctan(/3rfn/ adn),  cos 0dn  |odn|/7dn,  sin 0dn  (ov/n / 1 cxdn | )  /ij,, /yjn ,  (33) 

where  — 7r/2  <  6dn  <  tt/2.  Since  all  modes  have  the  same  c,  their  frequencies  are  proportional  to  the 
streamwise  wavenumbers, 


{lop,  los ,  udn}  =  { av ,  as,  \otdn\]  x  c  for  1  <  n. 


Thus, 


iod i  =  Los  —  ojp,  ujdN  =  Nujp  —  us,  udn  =  us  —  nojp  for  2  <  n  <  N  —  1. 


(34) 

(35) 
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From  (16),  (30)  and  (34),  one  may  show  that,  if  N  >  2, 

Ud(N-i)  <  up,  udN  <  up,  Up  <  udn  for  1  <  n  <  N  -  2.  (36) 

The  frequencies  of  the  first  N  —  2  difference  modes  are  higher  than  the  primary-wave  frequency.  But, 
those  of  the  last  two  difference  modes,  the  (N  —  l)th  and  Nth,  are  lower  than  the  primary  frequency. 

The  magnitude  of  the  unsealed  wave  elevation  is  O(e)  for  the  primary  wind  wave  and  0(e5)  for  the 
secondary  wave.  The  scalings  of  e  and  5  will  be  determined  later  (see  (69)  and  (70)).  The  difference¬ 
mode  elevations  become  0(ae5 )  in  water  since  we  only  consider  the  case  where  difference  modes  do 
not  satisfy  the  relation  (11). 

Equation  (6.5)  in  Lee  (2012)  shows  that  the  linear  growth  rate  of  a  wind  wave  is  of  0(a).  Thus, 
the  scaling  for  the  slow  time  t\  over  which  the  wave  growth  occurs  becomes 

ti  =  at.  (37) 

The  amplitude  functions  that  account  for  the  slow  growth  of  waves  are  denoted  by  Bt(ti)  for  i  = 
{p,s,dn}. 

Both  leading  and  second  order  velocity  potentials  (j)[  \z)  and  (j\  \z,t±),  for  i  =  {p,s,dl, . . .  ,dN}, 
satisfy  the  Laplace  equations.  The  solutions  for  the  primary  and  secondary  waves  become, 

=  cosh  [7 l(z  +  H )]  /  cosh(7iiL),  ^  =  0,  for  *  =  {p,  s},  (38) 

where  they  are  normalized  by 

c t>[  ^  +  a(j\  ^  =  1  at  2  =  0.  (39) 

For  the  difference  modes  the  solutions  are 

(i>dn  =  0,  4>dn  =  dn  coshfy dn(z  +  H)\Ddn(ti)/ cosh(7 dnH),  (40) 

for  1  <  n  <  N,  where 

dn  =  adnc2/[{G^fdn  tanh(7 dnH)  -  a2dn  (?}^dn  tanh(7 dnH)\,  (41) 

and  Ddn(t\ )  is  determined  by  nonlinear  critical-layer  interactions. 


Unsteady  flow  in  the  air 

Air  fluctuations  outside  the  critical  layer  are  also  governed  by  linear  equations  to  the  required  order  of 
approximation.  They  are  written  as,  in  the  limit  a  — )•  0, 
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The  magnitude  of  velocity  and  pressure  fluctuations  is  O(e)  for  the  primary  mode,  and  0(e5)  for  the 
secondary  and  difference  modes  in  the  air.  The  shape  functions  JJl(z,t\),  Vr(z,t\)  and  II^z, fi),  for 
i  =  {p,  s,  dl, . . . ,  dN},  are  related  to  <ht(z,  ti)  by  equations  (3.14)  and  (3.15)  in  Lee  (2012). 
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In  this  analysis,  the  nonlinear  jump  of  the  streamwise  fluctuation  velocity  across  the  critical  layer 
occurs  at  the  same  order  as  the  linear  jump.  Therefore,  the  shape  function  is  expressed  as,  for  i  = 
{p,s,dl,,..,dN}, 

Mz,t  i)  =  *ia\z)  +  ACti)^),  (43) 

where  will  be  determined  from  the  critical-layer  solutions,  and  4>(a^  and  <5^  satisfy  Rayleigh’s 

stability  equations 

(U  -  c)  (d2z  -  742)  4>’6)  -  U"$[a^  =  0  with  dz  =  d/dz,  (44) 

subject  to 

4>(“)  =  {i/'y1)dz(jil  ^  =  itanh(7 lH)  at  2;  =  0,  44“)  — >  0  as  z  —>■  00,  (45) 

4>(fe)  =0  at  2  =  0,  ->  0  as  2  — >  00.  (46) 

The  second  shape  function  is  the  second  linearly  independent  solution  of  Rayleigh’s  equation, 
which  can  be  obtained  by  the  method  of  variation  of  parameters  as  in  Reutov  (1980),  Drazin  &  Reid 
(2004)  and  Lee  (2012), 


= j,Ma)  -  < 


0,  2>2C 
$[a)  r&a))-2dt,  z<zc 


(47) 


where 

Jt  =  fZCm-*dz.  (48) 

Jo 

The  primary  mode  remains  linear  throughout  the  entire  analysis  and  Dp  =  0.  If  IV  >  2  or  up  <  ujs, 
the  secondary  mode  also  remains  linear  and  Ds  =  0.  All  difference  modes  are  assumed  not  to  satisfy 
the  free-surface-wave  dispersion  relation  (11),  therefore,  =  0. 

The  first  condition  at  2  =  0  in  (45)  was  obtained  from  the  kinematic  interface  condition  and  (38) 
with  (39),  and  could  be  viewed  as  normalization  conditions  of  4>p  and  4*5  for  the  linear  outer  solutions 
of  the  primary  and  secondary  modes.  From  (45)  and  (47),  we  obtain 


=  — i/  tanh(7jiL)  at  2  =  0, 


(49) 


for  1  =  {p,  s,di, . . . ,  dN}.  Note  that  4>j^  for  difference  modes  are  also  expressed  by  the  above  equation, 
which  can  be  considered  as  a  normalization  condition  for  4>^n  since  4>^  =  0. 

One  can  obtain  a  relation  that  expresses  the  modal  growth  rate  in  terms  of  4>jf4  and  Dt(t\).  For 
the  primary  and  secondary  modes,  1  =  {p,  s},  it  is  written  as 


2iy lB[/ (alcBl) 


tanh(y ZH)  + 


K  )\bi\2U"/U'c  +  A(*i)  =  Ri, 


(50) 


where 

/»oo 

Rt=  —  T  [|A^(a)|2  +  {%2  +  ^7(^-c)}|4>^|2]d2-7Jtanh(7Ji7)  +  (t/'/c)tanh2(7.8iL), 
Jo 


dz  =  lim 

e-5>0 


(51) 
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The  prime  denotes  the  differentiation  with  respect  to  the  relevant  argument,  i.e.  Bt  =  dBl/dti, 
U"  =  d2U/dz2,  and 

U'  =  U'(z  =  0),  (53) 

c  =  Uc  =  U(z  =  zc),  U'c  =  Uz(z  =  zc),  Ul(!  =  Uzz{z  =  zc),  (54) 

where  z  =  zc  is  the  location  of  the  critical  level.  The  bx  and  are  determined  by 

bt  =  ^a\z  =  zc)  =  lim  <&«(*),  C  -  b;  =  ( U'c/U'l/bt )  [dMa)\f-  (55) 

z^zf 

By  matching  the  outer  solution  with  the  critical-layer  solution  we  can  show  that  }\  —  bt  equals  i7r  as  in 
(86).  There  is  no  relation  corresponding  to  (50)  for  the  difference  modes  since  l1?  =  0  (see  Appendix 
B  in  Lee  2012).  In  Appendix  A,  it  is  found  that 

JA  =  lim  $J6>(z)  =  -i|JA|,  Jr  =  —iU'cRe(bi) / {'irU”\bl\2bl) ,  (56) 

Z->Zc 

for  i  =  {p,  s,d i, . . . ,  div},  where  Re  denotes  the  real  part  of  a  complex  variable.  For  numerical  evalua¬ 
tions  of  ./*,  the  above  formula  is  more  convenient  than  the  integral  form  in  (48). 


Outer  solution  expansion  about  the  critical  level 

The  velocity  and  pressure  field  near  the  common  critical  level  is  obtained  by  expanding  the  outer 
solutions  about  z  =  zr  as 


u  —  (c  +  U'cz  +  \U"z2  +  . . .)  = 


bx£lx  ( sin2  9t 

u" 

cos  9l  \  z 

U'c 

(ill  \z\—\  sin2  9t+bf  )  j+ 


[d.ja  cos  el) 

V  =  e6lBtblQl( sin0*)  [l/z  -  \  {U"/U'C)\  £x  +  c.c.  +  . . . , 

l 

w  =  -  eSli'ylBlbtQl  1  +  (U”/U'c)z  (in \z\  +bf  -1) 


£x+c.c.+  - . . , 


£%  +  c.c.  + 


p  —  Pp  +  Gz  =  yt  edxBxbAi( cos  9l)U'c  £x  +  c.c.  + 


(57) 

(58) 

(59) 

(60) 


where 

z  =  2;  -  zc,  (61) 

5p  —  1,  5S  —  ddn  —  d,  p  —  —  I;  £^dn  —  JdnBdm  Bp  —  Bs  —  0,  (02) 

for  1  <  n  <  N.  The  denotes  the  summation  of  all  2  +  N  modes,  i.e. 

X^(  )z  =  (  )p  +  (  )s  +  (  )di  +  •  •  •  +  (  )dN-  (63) 

The  last  term  in  braces  in  (57)  becomes  0  if  z  >  zc  and  Dx/(bx  cos  9X)  if  z  <  zc,  and  bx  becomes  6,+  and 

bt  above  and  below  the  critical  level.  The  above  solution  becomes  singular  at  z  =  zc,  therefore,  it  has 
to  be  rescaled  in  this  region. 

The  streamwise  velocity  jump  across  the  critical  layer  from  the  outer  solution  becomes, 

L.l'2'— +°° 


*Z=— oo 


secdt)  (U'J /U'JbAA  -bt)-Dt/bt 


£%  +  c.c.  +  . . . . 


(64) 


This  is  matched  with  the  jump  from  the  critical-layer  solution  to  derive  amplitude  equations. 
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CRITICAL-LAYER  EQUATIONS 


The  thickness  of  a  non-equilibrium  critical  layer,  where  the  growth  term  is  balanced  with  the 
linear  convection  effect,  is  of  the  same  order  as  the  small  growth  rate  as  shown  in  Reutov  (1980), 
Goldstein,  Durbin  Sz  Leib  (1987),  Goldstein  &  Hultgren  (1988),  Lee  (2012),  and  others.  Thus,  the 
vertical  coordinate  in  this  region  is  scaled  by 

C  =  (z  -  zc)/a  =  z/a.  (65) 

Viscous  terms  will  enter  the  critical-layer  momentum  equations  when  the  Reynolds  number  in  air  is 
scaled  by 

1  /Ra  =  o"3A,  (66) 

where  A  is  the  viscous  parameter  and  Ra  =  UaXa/ua  with  ua  denoting  the  kinematic  viscosity  of  air. 

The  continuity  and  momentum  equations  for  the  velocity  u  =  iu+jv+kw  and  pressure  p  expressed 
in  terms  of  the  scaled  variables  t\ .  X  and  (  become 


u 


■x  +  vy  +  w^/a  =  0, 


d  d  d  w  d  ~  d2 

’TWI+(U-C)dX+VTy+*Fr  X>\ 


{u,v,w}  =  -  \px,py,G+  -Pc 


(67) 

(68) 


The  magnitude  of  a  difference  mode  in  air,  that  is  generated  by  a  quadratic  interaction  between  a 
primary  wave  and  a  secondary  or  other  difference  mode,  becomes  of  the  same  order  as  the  secondary 
amplitude  when 

e  =  a3.  (69) 


Meanwhile,  the  nonlinear  effect  due  to  the  quadratic  coupling  of  the  secondary  and  lst-difference  modes, 
or  of  a  pair  of  difference  modes,  becomes  negligibly  small  when  5,  which  is  a  measure  of  the  amplitude 
ratio  of  the  secondary  to  primary  modes,  satisfies 


<y « i. 


(70) 


Introducing  (65)  and  (69)  into  the  outer  solutions  (57)-(60)  and  re-expanding  the  result  for  small 
(  show  that  the  critical-layer  solution  should  be  of  the  form 


f  u\ 


w 

\  p ) 


(  c  +  aU'c C  \ 
0 
0 

\  Po-Gz  ) 


+  cR 


/  r/1)  \ 

r/1) 


+  cG  In  a 


(  \ 

w(2  L) 

awW 
(2  L) 


+  a° 


l  ■yl2)  \ 


aw ^ 


\  crpW  J 


+ 


(71) 


y  ap W  J  y  apyAiJ)  J 

It  is  convenient  to  introduce  the  critical- layer  amplitude  At, 

A  =  blQtBl  for  i  =  {p,  s ,  dN},  Adn  =  -b*dnn*dnB*dn  for  1  <  n  <  N- 1, 


(72) 


where  Bt  is  the  amplitude  in  the  outer  solutions,  and  D*  is  defined  by  (62).  Then,  the  leading-order 
term  that  matches  onto  the  outer  solution  (57)-(60)  is  written  as 


( u^-\ufe 


V 


2  c  ' 
(1) 


=  1 


( 


fip/otp 

-1 


Qp  £p-\-5i 


Qdn^dn+C.C.,  (73) 
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Adn£dn  “1“  C.C..  (74) 


( 


P 


(1) 


-17 V 
Ucoip/ 7 p 


ApSp  +  5 


-17 5 

U'cOtshs 


N 

as£s + 5  yy 

n=  1 


17  dn 

UcOtdn/  7 dr 


The  wavenumbers  are  defined  by  (12),  (30)  and  (32).  The  and  £s  are  given  in  (29),  and  £dn  is 
defined  as 

£dn  =  exp[i (adnX  +  f3dny )]  for  1  <  n,  (75) 

that  becomes  £dn  =  £dn  it  1  <  n  <  N  —  1,  and  £dN  =  £d./v  with  £dn  given  in  (29).  The  second-order 
critical-layer  solution  must  include  all  harmonics  generated  by  the  nonlinear  terms  in  the  momentum 
equations.  They  become,  up  to  0(5), 


(2)  _  A(2)  C  ,  /7(2)  J  yd2)  c2  ,  r 

q[  —  Qp  £p  +  Q o  +  Q2P  £p  +  5 


N+ 1 


n=l 


+  C.C., 


(76) 


where  and 


denote 


Qf  =  [uf\vf\ wf’.pf}  for  7={p,0,  2j>,  s,j>  +  s,  <il,  . ...  d(jv+l)>.  (78) 

By  substituting  (71),  (73),  (74)  and  (76)-(78)  into  the  continuity  and  momentum  equations,  a  system 
of  critical-layer  equations  is  derived  for  q[ 1‘>  and  as  given  in  Appendix  B. 
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AMPLITUDE  EQUATIONS 


Matching  of  critical-layer  jumps 

The  critical-layer  equations  in  Appendix  B  are  analytically  solved  to  obtain  the  velocity  jumps  across 
the  critical  layer  as  in  Goldstein  &  Choi  (1989),  Goldstein  &  Lee  (1992),  Wu  (1992),  Wu,  Lee  &  Cowley 
(1993),  Lee  (2012),  and  others.  The  following  normalized  variables  are  introduced, 

at  =  at/a,  fa  =  fa/ fa  fa  = 'y t/fa  (79) 

t  =  k(aU'ct\  —  t0) ,  (  =  (/&  —  ( o,  X  =  X  —  XQ,  A  =  X/(ak3Urc),  (80) 

At  =  (a/ M  /U’c/k3)Al ei«d^o+Cot))  jj( 2)  =  {a2  M  /  (f32  k3)}  uf  ,  (81) 

where  the  normalization  parameters  a,  fa  k  and  M  and  the  coordinate  origin  shift  (Q  are  chosen  in 
Appendix  C,  and  tQ  and  Xa  will  be  chosen  later. 

When  N  >  2,  the  velocity  jumps  of  the  primary,  secondary  and  N  difference  modes  become 

/OO 

ujfd(  =  m{a./P)\jM/U'c  U'fafa/afa  Afat)  for  j  =  {p,  s},  (82) 

-OO 


ufldQ  = 


Ki  +  , 

A-n  +  Ffa  +  G+ , 

Ait  +  G„  , 


1  fan  <N- 2, 
n  =  N—  1, 
n  =  N, 


where  A±,  Fn  and  are  given  in  Appendix  D. 

The  difference  mode  jumps  (83)  are  rewritten  as,  for  example,  if  N  =  2, 


(83) 


(t\Adi)  +  F1  ( t\Ap ,  A*)  +  G/  ( t\A *,  A^), 
A-2~ (J\Ad2)  +  G2  (t\Ap,  A^i), 


n  =  1, 

n  =  2. 


(84) 


Two  difference  modes  are  generated  by  a  2  +  2  mode  interaction  as  shown  in  Figures  2(e,  f)-4(e,  f). 
When  N  =  3,  equation  (83)  for  a  2  +  3  mode  interaction  becomes 


A-i  (t\A(ii)  +  Fj  (t\Ap,  As), 

<  A2  (t\Ad2)  +  F2  (t\Ap,  Adi)  +  G2  it\A*p,  Adz), 
A3" (t\ Adz)  +  G3  (t\Ap,  Ad2), 


n  =  1, 
n  =  2, 
n  =  3. 


(85) 


Diagrams  of  2  +  3  mode  interactions  are  illustrated  in  Figures  2(g,  h)-4(g,  h). 

Matching  between  the  critical-layer  solutions  (71)  and  the  outer  solutions  (57)-(60)  requires  that 
the  velocity  jumps  (82)  and  (83)  from  the  critical-layer  solutions  at  0(cr35l)  must  be  equal  to  the 
corresponding  velocity  jumps  (64)  from  the  outer  solution.  By  matching  them,  we  obtain 


=  ivr, 


(86) 
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DP  =  DS  =  0,  DdnBdn  =  —  /3bdnadn/ (cry dn\JU'cM)  Fdn(t)  for  1  <  n  <  N,  (87) 

where  Fdn  denotes  the  nonlinear  jumps  without  the  terms  on  the  right  sides  of  (83).  The  normalized 
outer  amplitude  Bt  is  related  to  the  normalized  critical-layer  amplitude  A%  by 

A,,  =  bintBl  with  Bt  =  (\/ Mjui/k3)Blei&'{Xo+ <«>*>.  (88) 

The  above  equations  (86)— (88)  with  the  relation  (50)  completely  determine  the  evolution  of  amplitudes. 
They  can  be  combined  to  derive  more  explicit  amplitude  equations  as  will  be  summarized  below. 


Primary  and  secondary  mode  amplitudes 

The  primary  mode,  whose  magnitude  of  0(< r3)  is  larger  than  the  other  modes,  remains  linear  throughout 
the  whole  stage.  When  TV  >  2,  or  ujp  <  us.  the  secondary  mode  of  magnitude  of  0(a36)  also  remains 
linear.  From  the  relation  (50)  and  the  matching  conditions  (86)  and  (87),  we  can  show  that 

Ap  =  apexp  (Kpt),  As  =  asexp(nst).  (89) 


The  linear  growth  rates  np  and  ns  become 

7rc£/'/dl|&j|2 

- - - - - l _ l _  —  YCX 

2k(3U'(?A/l  tanh(7*77) 


2k/3U'cA(l  tanh(7,i7) 


-Co  for  z  =  {p,s}, 


with  /3,  k,  Co  and  Rt  in  (C-l)-(C-3)  and  (51).  The  initial  amplitudes  ap  and  as  are 

at  =  (A[0^/Al °o)ex.p  («j  -  l)kt0  +  i(dj  -  d30)(Xo  -  kt0C0)  for  i  =  {p,  s}, 

where  A[°^  is  a  non- normalized  initial  amplitude  of  Au 

A  — 

S±%  —  , 

and  Kt  is  the  non-normalized  linear  growth  rate, 

fq  =  (dAl/dti)/Al  =  akU'c{Kr  -  id,C0). 

The  coordinate  origin  shifts  t0  and  XQ  introduced  in  (80)  are  chosen  to  satisfy 

(y/Mjul/k3)A.^exp  ktQ  +  ia30(Xo  -  kt0(o)  =1. 


By  putting  (C-2)  and  (C-3)  into  (90),  and  from  (91), 

k30  =  a30  =  L  (95) 

In  Lee  (2012),  we  showed  that  kp  and  ap  could  be  normalized  to  unity  by  selecting  the  properties  of 
the  primary  and  secondary  waves  of  each  set  of  interactions  as  normalization  parameters.  The  as  could 
also  be  normalized  to  be  one.  However,  in  this  paper,  the  general  forms  in  (89),  with  normalization 
parameters  chosen  from  the  properties  of  a  wind  wave  with  30°  propagation  angle  will  be  kept  for 
straightforward  comparisons  between  numerical  results  of  various  sets  of  interactions. 

When  the  secondary- wave  frequency  is  lower  than  the  primary  frequency,  or  N  =  1 ,  the  2  +  1  mode 
interaction  produces  a  nonlinear  effect  on  the  secondary  wind  wave.  Its  amplitude  is  then  governed  by 
an  integro-differential  equation. 
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Amplitude  equations  of  difference  modes  when  A  =  0(1) 

Equation  (87)  shows  that  the  amplitudes  of  the  first  N  —  2  difference  modes  grow  exponentially, 

A dn  =  adnexp(Kdnt)  for  1  <  n  <  N- 2.  (96) 

Their  growth  rates,  ndn,  and  initial  amplitudes,  adn,  become 

Kdn  =  K*s+  nKp,  (97) 

®dn  P dn^'P®'  d(n-l)^ d  with  ddQ  ~0'si  (98) 

where  ipdn  and  Id^  are  given  by 


AdrJ^dn  _  .  b, dn  T30 ®-dn  I'^dnll^dnl 

«30 Idn  l^3oll^3oP  \bdn\  O30 Idn  I^3oII^3oP 


(99) 


rW 

Ld 


Vnld(n-X)  f0°°  4-}(r)  r2  dr 

-C^J  f0°°  fo°°  4-Vi,  r2)  e-^P+<)(Tl+r2)  Tl  dTldT2 


(100) 


with  z>,  un,  Cc-\  c/A ,  E i”'1  and  E^}  in  (D-4)-(D-7).  The  relation  (56)  was  used  to  obtain  (99). 

However,  the  amplitudes  of  the  last  two  difference  modes,  the  ( N  —  l)th  and  IVth,  are  determined 
by  coupled  integral  equations 


A-d(N-i)  —  ad(jv-i)exP(Kd(jv-i)  f)  +  AdN), 


(101) 


A<In  —  ^PdN^N {t\Ap,  H^cjv-i)),  (102) 

with  ud(N_V)  defined  by  (98)  and  G±  given  by  (D-3)  in  Appendix  D.  The  above  equations  determine 
Ad(N_ i)  and  AdN  at  a  given  i. 

From  (101)  and  (102),  one  may  obtain  that 


Ad(N-i)  ->  arf{Ar_1)exp(Kd{Ar_1)t),  AdN  -»  adNexp(KdNt),  as  t  -»  -oo,  (103) 

where  K^jv-i),  Krfjv  aRd  arf(iv-i)  are  defined  by  (97)  and  (98),  and  adN  is  equal  to  ro  in  (125). 

Equations  (36)  and  (96)  show  that  the  first  N  —  2  difference  modes,  whose  frequencies  are  higher 
than  the  primary-wave  frequency,  grow  exponentially.  The  growth  rate  of  the  nth-difference  mode  is 
increased  by  kp  and  its  frequency  is  decreased  by  up  compared  to  those  of  the  preceding  (n  —  In¬ 
difference  mode.  As  shown  in  (103)  and  (115),  the  growth  rates  of  the  ( N  —  l)th  and  IVth  difference 
modes  in  the  initial  exponential-growth  stage  are  also  larger  by  kp  than  the  ( N  —  2)th  and  ( N  —  l)th 
modes,  respectively.  However,  the  growth  rates  of  these  last  two  difference  modes  can  become  much 
greater  than  their  initial  growth  rates  when  i  becomes  large.  The  frequencies  of  the  ( N  —  l)th  and 
IVth  difference  modes  are  lower  than  the  primary  frequency  as  given  in  (36). 

The  amplitudes  of  all  2  +  N  modes  ( N  >  2)  are  now  determined  by  (89),  (96),  (101)  and  (102). 
It  is  noted  that  only  initial  amplitudes  of  the  wind  waves,  ap  and  as,  can  be  freely  chosen  as  in  (89). 
All  initial  amplitudes  of  nonlinearly  generated  difference  modes,  adn ,  are  determined  by  (98).  The 
governing  equations  become  much  simpler  in  the  inviscid  and  large-viscosity  limits  as  will  be  given 
later.  We  will  also  obtain  the  series  solutions  of  the  ( N  —  l)th  and  IVth  difference-mode  amplitude 
equations. 
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Outer  amplitudes  Bp,  Bs  and  Bdn 

From  (62),  (87)  and  (88),  the  primary  and  secondary  amplitudes  in  the  outer  solutions,  Bp  and  Bs, 
that  are  normalized  by  the  second  equation  in  (88),  become 

Bp  =  Ap/bp,  Bs  =  As/bs,  (104) 

where  Ap  and  As  are  given  by  (89).  The  difference  mode  is  solely  generated  by  the  critical-layer 
interaction,  and  Bdn  always  appears  multiplied  by  Ddn  in  the  outer  solution  (57)-(60).  From  (88)  with 
(62),  the  term  DdnBdn  becomes 

BdnBdn  —  Adn/ {bdnJdn)  ■  (105) 

The  outer  solutions  in  the  air,  the  perturbations  in  the  water  and  the  wave  elevation  are  now  completely 
determined. 
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INVISCID  AMPLITUDE  EQUATIONS 


In  the  inviscid  limit  where  A  =  0,  (96)  becomes 


Adn  =  a^u) exp(Kdnt)  for  1  <  n  <  N- 2, 

(106) 

with 

(inv)  *  (inv)  T(n) 

adn  ~  VdnaPa din-1)1  dinv’ 

(107) 

41  =  44  =  0)  =  C<”>  /  (**,  +  k;)3  . 

(108) 

Equations  (101)  and  (102)  are  simplified  as 

/»oo 

Ad(N_1}  =  aJ(X)1)exp(Kd(JV_1)f)  +  e2(Reft^  /  (r)  (*  -  r)dr, 

J  0 

(109) 

eKpt  /  e^(1_l/jv)Kpr  t2  Ad(JV_i)  (t  -  r)dr, 

Jo 

(110) 

where 

K%nv)  =  e~{1~UN)KpT  [ r 2  -  3r/(i/JVRe/cp)  +  3/(uNReKp)2] 

_e-(KP+vNK*P)T  ^r2  +  3T /(z/^ReAvp)  +  3/(i/iVRe/ip)2]  , 

(111) 

*7* Rnv)  =  _l|ap| 2^JV_1)C,ff“l)  (^Tp  -  |C,£V_1)) 

(112) 

=  -aP<PdNC{aN)  (uN%  -  \C^N)S)  . 

(113) 

From  (D-4)  and  (16), 

uN  =  N  —  v  =  N  —  as/ap  and  0  <  uN  <  1. 

(114) 

Equations  (89),  (106),  (109)  and  (110)  determine  the  inviscid  amplitudes.  The  series  solutions  of  (109) 
and  (110)  will  be  given  later.  We  will  also  find  the  conditions  when  the  coefficients  tpanv ^  and  (pjnv^ 
vanish. 

The  amplitudes  A^iv-p  and  AdN  in  (109)  and  (110)  become 

Ad{N-i)  ->  aJ{7)1)exp(Kd(JV_1)f),  AdN  ->  a{™v)ex.p(KdNt),  as  t  ->  -oo,  (115) 

where  }  is  given  by  (107),  and  ad™^  equals  ro  defined  later  by  (125). 
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SOLUTIONS  OF  VISCOUS-LIMIT  AMPLITUDE  EQUATIONS 


In  the  large-viscosity  limit  where  A  — >  oo,  (96)  can  be  written  as 

Adn  =  a*fm)exp(ftdnt)  for  1  <  n  <  N- 2, 


where 


{vlim )  _  *  ( vlirn )  An) 

dn  r  dn  d(n— 1)  dviim ’ 


(116) 

(117) 


An)  =  Um4n)=  Cin) 

dviim  A-S-OO  “ 


7d(n_i)/{Aa^n(l-i/n)}-  /“JEi")(ri,T2)ridTi(ir2 


(118) 


The  viscous-limit  amplitude  equations  for  the  (IV  —  l)th  and  Arth  difference  modes  are  obtained  by 
taking  the  A  — >  oo  limits  of  (101)  and  (102)  using  similar  methods  in  Wu  et  al.  (1993)  and  Lee  (2012). 
These  equations  are  solved  to  obtain  the  following  analytic  solutions 


id(JV-l) 


=  o£SexP(«d(jv-ii*)/(1  -  6 


^d(iV-l) 


'd(AT-l) 1 


_  Avlim) 


;2(R eKp)t^ 


A 


d(N-i) ) 


(119) 


where 


Avlim) 

ra 


CLp 


*<P. 


d(N-i) 


~  (vlim)  -=,+ 

<Pb 


/v  (vlim)  Pv— 

T’ft  —  ~apAdN^N1 


(120) 


=  T  C*a± (ra)  7P/ {AdpZ/„(l  ±  r'n)2}  -  cf^  /0°°  ^i±j(ri,  r2)  n  dTidr2 


An)i 


(121) 


By  renormalizing  the  amplitudes  and  time,  we  may  find  that  A  is  not  an  independent  parameter  for 
the  above  solutions.  In  the  large-viscosity  limit,  the  evolution  of  all  amplitudes  is  described  by  simple 
algebraic  equations  (89),  (116)  and  (119). 
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SERIES  SOLUTIONS  OF  AMPLITUDE  EQUATIONS 


Both  inviscid  and  finite-viscosity  amplitude  equations,  in  (101),  (102),  (109)  and  (110),  possess  the 
following  absolutely-convergent  series  solutions, 

OO  OO 

Ad(N-i)  =  ^  sn  exp[(2raRe/%,  +  Kd(JV_1})  t],  AdN  =  ^  rn  exp[(2raRe/tp  +  ndN )  t],  (122) 

n=0  n= 0 


where  Kdn  is  defined  by  (97).  The  recursion  relations  are  obtained  as 


(inv)  -r  \  n  •  r 

s0  —  ad(N-l )  A  —  0,  So  —  ad(N- 1)  n 

A  /  0, 

(123) 

Sn  =  (1/4 )<P%nv)vfr(n  -  1  +  £)_3(n  -  1  +  £  +  i/iV)_3(ReKp)“ 

3fl,n- 1  sn_i  for  n  >  1, 

(124) 

rn  =  (1/4 )<p£nv\n  +  C)_3(Re^p)-3  K  sn  for 

n  >  0, 

(125) 

where 

£  =  +  k*)/  (2R  enp). 

(126) 

In  the  inviscid  limit,  gn-i  and  hn  become 

9n- 1  =  hn  =  1  if  A  =  0, 

(127) 

and  they  are  given  by  (E-l)  and  (E-2)  in  Appendix  E  for  the  finite-viscosity  case. 

The  recursion  relation  (124)  is  summed  to  obtain 

=  (#",)/4)n  k/rsk,,)3”  [©„« + vN)„r3  (n:Ui  sm.d «».  <128) 

where  (£)n  =  T(£  +  n)/T(£)  denotes  the  generalized  factorial  function,  and  T(£)  is  the  gamma  function. 

Following  the  method  used  in  Wundrow,  Hultgren  &  Goldstein  (1994),  we  can  show  that  the 
asymptotic  behavior  of  (122)  when  t  — >  oo  becomes 


Sasyex  p 


x  exp 


6 


{ 

{ 


Kd(N-i)  -  (2£  +  VN  -  7/6) Re/ip}  t 
(I/jV/Ite«p)3^ont’)/4}  1  exp{(Re/tp/3)f}  , 


(129) 


AdN  (t)  ~  rasyex p  [i(Im/ip)t]  Ad(Ar_1}  (t) ,  (130) 

where  s  and  rasy  are  given  in  (E-7)  and  (E-8),  and  Im  denotes  the  imaginary  part  of  a  complex 
variable. 

Similar  to  previous  analyses  by  Goldstein  &  Lee  (1992),  Wundrow  et  al.  (1994)  and  Lee  (2012), 
the  amplitudes  Ad(jv-i)  and  AdN  grow  like  the  exponential-of-an-exponential  as  t  — >  oo.  It  is  worth 
noting  that  the  viscous  effect  only  appears  in  the  coefficients  sasy  and  rasy  in  (129)  and  (130).  Their 
f-dependent  parts,  or  Ad(N_V)/ sasy  and  AdN/rasy,  are  independent  of  A. 
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SPECIAL  SOLUTIONS  OF  AMPLITUDE  EQUATIONS 

Generation  of  plane  second-difference  modes  by  2  +  2  mode  interactions 

By  putting  [3dn  =  0  and  using  (11),  (12)  and  (30),  we  may  show  that  a  2  +  2  mode  interaction  can 
generate  a  second-difference  mode  whose  motion  is  two-dimensional,  or  0d2  =  0,  if  ap  and  as  satisfy 

a2  =  ^1  +  yj  1  +  16 ms(mpa^  -  1  )a^j  /(2 ms),  (131) 

where  m,  =  c4/[G2  tanh2(7ji7)]  for  i  =  {p,  s}.  The  combinations  of  (0P,  0S)  that  satisfy  (131)  are 
plotted  as  the  6d2  =  0  contour  line  in  Figure  6(b)  and  also  as  red  curves  in  Figures  7(a,  b). 

One  can  further  show  that  the  plane  2nd-difference  mode  is  a  wind  wave,  if 

(4  -  mv/ms)f2Vp  -  8fVp  +  6v2  +  (1  —  2 /f2)rvp  +  (1/f2  -  l)/4  =  0,  (132) 

where  vp  =  ^/rrisap.  r  =  tanh(7sFT)/ tanh(a2Dl7)  and  07 D  =  (G/c2)  tanh(a2D-ft).  In  the  deep  water 
limit,  where  tanh^H)  — >  1  and  tanh(a2D-ff)  — >  1,  (132)  becomes 

ap  =  (5  +  \/l3)G/(6c2).  (133) 

With  parameter  values  in  (14)  and  (134),  the  condition  is  satisfied  by  the  primary  and  secondary  waves 
with  0p  «  45.8°  and  0S  «  57.6°.  The  above  results  (131)-(133)  are  independent  of  A. 

If  the  plane  2nd-difference  mode  is  not  a  wind  wave,  the  amplitudes  of  both  difference  modes, 
Adi  and  ^4^2 1  grow  exponentially  and  they  are  determined  by  (103)  or  (115)  with  N  =  2.  When  the 
condition  (132)  or  (133)  is  satisfied,  Ad2  will  be  determined  by  an  equation  that  is  similar  to  (5.12)  in 
Lee  (2012). 

When  as  =  Nap 

If  as  =  Nap,  the  critical-layer  interaction  generates  N  difference  modes  with  wavenumbers  given  by 
(30).  Since  the  streamwise  wavenumber  of  the  IVth-difference  mode  is  zero,  its  outer  solution  is  not 
singular  and,  as  a  result,  the  difference-mode  or  nonlinear-effect  generation  mechanism  is  stopped.  If 
N  =  1  and  f3s  =  f3p,  the  IVth-difference  mode  becomes  trivial  for  the  primary  oblique  wave  interacting 
with  itself  (Goldstein  &  Leib  1989).  In  other  cases,  i.e.  if  IV  =  1  and  f3s  =  —  /3P,  or  if  N  >  2,  the  IVth- 
difference  mode  becomes  steady  and  spanwise  periodic  (Goldstein  &  Choi  1989).  The  exponentially- 
growing  amplitudes  of  the  first  IV  —  1  difference  modes  are  determined  by  (96)  and  (103).  When  N  =  0, 
the  plane  secondary  wave  grows  exponentially,  but  the  difference  mode  is  not  generated  as  given  in 
Appendix  F.2  in  Lee  (2012). 

When  ipa'nv>  =  0  and  <^nv>  =  0 

The  coefficients  (pa and  1  in  the  inviscid  amplitude  equations  (109)  and  (110)  can  become  zero. 
In  Figure  7(a),  the  combinations  of  ( 0P ,  0S)  for  =  0  are  plotted.  The  contour  lines  of  tp<£nv^  =  0 

are  plotted  in  Figure  7(b)  for  the  parameter  values  in  (14)  and  (134).  Border  lines  for  various  2  +  N 
mode  interactions,  which  were  shown  in  Figure  5,  are  replotted  for  1  <  N  <  7  as  green  dotted  curves. 
The  red  solid  curve  depicts  the  (0P,  9S)  combinations  that  generate  plane  second-difference  modes.  It 
is  the  same  curve  as  the  0d2  =  0  contour  line  in  Figure  6(b).  Because  of  (112),  po'nv'1  =  0  if  ip<£nv^  =  0. 
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Figure  7.  Contour  lines  of  (a)  (fa™'*  =  0  and  (b)  p^nv'1  =  0  in  the  (9p,  6S)  plane  plotted  as  solid  curves. 
Border  lines  for  various  2  +  N  mode  interactions,  shown  in  Figure  5,  are  replotted  for  1  <  N  <  7  as 
green  dotted  curves.  Red  curves  are  identical  to  the  0d2  =  0  contour  line  in  Figure  6(b).  The  straight 
blue  line  in  region  (i&)  in  (a)  satisfies  the  condition  (F.l)  in  Lee  (2012). 


Thus,  red  and  black  solid  curves  are  identical  in  both  Figures  7(a,  b).  Blue  solid  curves  depict  the 
conditions  when  nu'1  =  0  but  (p ^nv'1  ^  0. 

If  =  0  and,  thus,  ipa'nv'>  =  0,  the  amplitude  grows  exponentially  as  given  by  the  first 

equation  in  (115),  but  Ac]iN  becomes  trivial.  If  ipa™^  =  0  but  cp^'nv^  /  0,  both  A^n-i)  and  A(jn  grow 
exponentially  as  determined  by  (115).  On  the  straight  line  in  region  (4)  in  Figure  7(a),  both  secondary 
and  difference  modes  of  a  2  +  1  mode  interaction  grow  exponentially  as  shown  in  Appendix  F.l  in  Lee 
(2012). 
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NUMERICAL  RESULTS 


All  numerical  results  are  obtained  by  using  the  logarithmic  profile  (7)  and  the  parameter  values 
listed  in  (14)  and 


a  =  0.0012. 


(134) 


The  quantities  are  non-dimensionalized  by  Ua  and  Aa  given  in  (8).  With  (4)  and  (14),  the  Charnock’s 
(1955)  constant  becomes  aCH  =  gz0/v%  =  G/ 9.025  =  5.54  x  10-3.  The  parameter  values  are  chosen  to 
show  the  characteristic  behaviors  of  the  solutions  of  amplitude  equations. 

Once  the  wave  propagation  angles  of  the  primary  and  secondary  wind  waves,  9P  and  9S ,  are  selected, 
the  wavenumbers,  frequencies  and  other  propagation  angles  are  determined  by  (11)— (13),  (30)  and  (32)- 
(34).  The  number  of  difference  modes  generated  by  a  2  +  N  mode  interaction  is  determined  by  (15). 
The  viscous  parameter  A,  defined  by  (66)  and  (80),  is  treated  as  a  free  parameter.  The  linear  shape 
functions  4>pl\  44^  and  4>^  are  obtained  by  numerically  solving  (44).  This  solution  is  used  to  evaluate 
Rt,  bt  and  Jt  by  (51),  (55)  and  (56).  All  parameters  appearing  in  the  amplitude  equations  are  now 
determined. 

The  amplitude  equations  are  solved  subject  to  initial  conditions.  As  mentioned  before,  only  initial 
amplitudes  of  wind  waves  can  be  arbitrarily  chosen.  Those  of  the  primary  and  secondary  wind  waves 
are  chosen  as 

ap  =  1,  as  =  0.001,  (135) 


for  all  numerical  results  in  this  section.  The  initial  amplitudes  of  all  difference  modes  in  a  2  +  N  mode 
interaction  are  determined  by  (98),  (107)  or  (117)  for  the  finite-viscosity,  inviscid  or  viscous-limit  case, 
respectively,  if  N  >  2. 

When  N  >  2,  the  primary  and  secondary  wind  waves  remain  linear  as  in  (89).  The  first  N  —  2 
difference-mode  amplitudes,  that  also  grow  exponentially,  are  determined  by  (96),  (106)  and  (116).  The 
series  solutions  given  in  (122)  are  used  to  compute  the  amplitudes  of  the  ( N  —  l)th  and  Nth  difference 
modes  in  the  finite-viscosity  and  inviscid  cases.  For  the  results  presented  in  this  section,  summations 
of  less  than  150  terms  were  sufficient  to  obtain  convergent  solutions.  The  integral  equations  (101), 
(102),  (109)  and  (110)  can  also  be  numerically  solved  by  using  the  procedure  developed  in  Goldstein 
&  Lee  (1992)  and  Lee  (1997,  2012).  In  the  viscous  limit,  the  last  two  difference-mode  amplitudes  are 
computed  from  the  closed- form  solutions  (119). 

Figure  8  illustrates  the  evolution  of  amplitudes  in  an  inviscid  2  +  2  mode  interaction.  The  propa¬ 
gation  angles  of  the  primary  and  secondary  wind  waves  are  chosen  to  be  9P  =  30°  and  9S  =  —57°  as 
denoted  by  the  letters  P  and  S  in  Figures  2(e,  f).  The  nonlinear  interaction  generates  two  difference 
modes,  denoted  by  D\  and  D2,  whose  propagation  angles  are  9 41  ~  —79°  and  64 2  «  83.5°.  The  wave 
system  in  Figure  1(b)  corresponds  to  this  case. 

The  linear  primary  and  secondary  amplitudes,  Ap  and  As,  are  shown  in  Figure  8(a).  The  amplitudes 
Adi  and  Ad2  of  the  first  and  second  difference  modes,  computed  using  the  finite- viscosity  series  solutions 
(122),  are  plotted  in  Figure  8(b)  as  solid  curves.  The  initial  exponential  solutions  (115),  which  are  equal 
to  the  first  terms  in  the  series  solutions,  are  shown  as  dot-dashed  curves.  The  asymptotic  solutions 
(129)  and  (130)  are  plotted  as  dotted  curves.  The  amplitudes  Adi  and  A 42  grow  exponentially  in  the 
linear  stage  where  t  <  0.5  although  they  are  nonlinearly  generated.  Both  difference  modes  grow  faster 
than  the  primary  and  secondary  wind  waves,  and  grows  faster  than  Adi  hi  the  initial  linear  stage. 
Once  the  primary  amplitude  reaches  a  certain  level,  say,  at  t  ~  0.5,  it  starts  to  affect  the  growth  of 
difference  modes.  In  the  later  part  of  the  parametric-growth  stage  where  t  >  2.5,  Adi  and  Ad2  grow 
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Figure  8.  Inviscid  2  +  2  mode  critical-layer  interaction  when  0P  =  30°,  6S  =  —57°  and  A  =  0.  (a) 
In  \Ap\  (black),  ln|As|  (green);  (b)  ln|v4di|  (blue),  In  |Ad2|_(red);  (c)  Re(A'dl/Adi)  (blue),  Re(A'd2/Ad2) 
(red);  (d)  lm.(A'dl/ Adi)  (blue),  lm.(A'd2/ Ad2)  (red)  versus  t.  Solid  curves,  0(l)-viscosity  solutions;  dot- 
dashed  curves,  initial  exponentially-growing  solutions;  dotted  curves,  asymptotic  solutions.  Blue  and 
red  dotted  curves  are  indistinguishable  and  only  red  dotted  curves  are  visible  in  (c,  d). 
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almost  identically.  They  also  grow  faster  than  exponentially.  In  fact,  they  grow  like  the  exponential-of- 
an-exponential  as  solid  curves  of  the  finite-viscosity  solutions  are  indistinguishable  from  dotted  curves 
of  the  asymptotic  solutions  in  Figure  8(b).  Transition  from  the  initial  exponential- growth  stage  to  the 
later  exponential-of-an-exponential-growth  stage  is  shown  to  be  gradual. 

The  real  and  imaginary  parts  of  A'dn/Adni  plotted  in  Figures  8(c,  d),  represent  the  growth  rate 
and  wavelength  reduction,  respectively.  The  real  part  also  represents  the  energy  transfer  rate  since 
(d\A,in\2 /dt) /\Adn\2  =  2Re(A'dn/ Adn).  The  0(l)-viscosity  solutions  are  shown  as  solid  curves  and 
the  asymptotic  solutions  are  plotted  as  dotted  curves.  In  the  figures,  blue  dotted  curves  for  the  1st- 
difference  mode  are  not  distinguishable  from  red  dotted  ones  for  the  2nd-difference  mode.  The  former 
is  invisible  since  it  is  completely  covered  by  the  latter.  In  the  initial  exponential-growth  stage,  the 
1st  and  2nd  difference  modes  have  constant  Re(Adn/Adn)  values  of  2.64  and  3.64.  In  the  entire  stage, 
Re(A!dl/ Adi  —  Kdl)  and  Re(A'd2/Ad2  —  nd2)  are  positive.  As  a  result,  |An|  and  |Ae|  grow  monotonically. 
The  constant  values  of  1  and  1.64  for  R e(A'p/Ap)  and  Re(A/S(Ms)  are  not  plotted  in  the  figure. 

Results  of  the  inviscid  2  +  2  mode  interaction  when  9P  =  30°  and  9S  =  50°  are  shown  in  Figure  9. 
The  difference- mode  propagation  angles  become  9di  ~  71.3°  and  9d2  ~  —34.6°.  Four  modes  involved 
in  this  interaction  are  denoted  by  P,  S,  D\  and  D 2  in  Figures  2(e,  f).  The  nonlinear  growth  effects 
start  to  become  large  near  t  ~  2.  In  the  initial  linear  stage  where  Adi  and  Ad2  grow  exponentially,  the 
signs  of  Re(A/dl/Arf1  —  ndl)  and  Re(Arf2/ A^2  —  remain  unchanged.  Since  they  are  negative  in  this 
case,  the  initial  nonlinear  effects  tend  to  reduce  the  growth  rates  as  one  may  observe  in  Figure  9(c). 
However,  the  continuing  nonlinear  effect  enables  the  reverse  of  their  signs  for  the  eventual  exponential- 
of-an-exponential  growth  of  amplitudes.  The  finite-viscosity  solutions  shown  as  solid  curves  are  not 
distinguishable  from  the  asymptotic  solutions  plotted  as  dotted  curves  when  t  >  5. 

Effects  of  the  viscous  parameter  A  on  the  growth  of  the  first  and  second  difference  modes  are  shown 
in  Figures  10(a,  b)  for  the  same  values  of  9p  =  30°  and  9S  =  —57°  as  in  the  inviscid  case  given  in  Figure 
8.  The  value  of  A  is  varied  from  0  (inviscid)  to  10°.  The  finite-viscosity  solutions  are  plotted  as  solid 
curves  along  with  the  corresponding  viscous-limit  solutions  (119)  as  dashed  curves.  The  asymptotic 
solutions  computed  by  (129)  and  (130)  are  plotted  as  dotted  curves.  The  linear  primary  and  secondary 
wave  amplitudes,  shown  in  Figure  8(a),  are  not  affected  by  A.  The  finite-viscosity  solutions  when 
A  =  10  plotted  as  red  curves  are  almost  identical  to  the  inviscid  ones  plotted  as  black  curves  in  Figures 
10(a,  b).  It  is  shown  that  the  onset  of  the  nonlinear  growth  is  delayed  as  A  is  increased,  this  can  be 
identified  as  the  departing  point  from  the  initial  straight  lines  (not  plotted  here,  but  similar  to  the 
dot-dashed  lines,  for  example,  in  Figure  8(b)).  We  also  note  that  Adi  and  A([2  grow  almost  identically 
in  the  nonlinear  growth  stage  for  all  A  values  considered.  The  inviscid  and  all  finite-viscosity  solutions 
approach  the  exponential-of-an-exponential  asymptotic  solutions  as  t  — )•  00.  Even  for  large  values  of 
A,  the  evolution  of  Adi  and  Ad2  in  later  nonlinear  growth  stage  is  approximately  described  by  the 
asymptotic  solutions  plotted  as  dotted  curves.  It  is  worth  pointing  out  that  there  is  no  arbitrary 
parameter  in  the  asymptotic  solutions  (129)  and  (130). 

As  A  becomes  large,  the  0(l)-viscosity  solution  plotted  as  solid  curves  in  Figure  10  approaches 
the  viscous- limit  solution  plotted  as  dashed  curves,  especially,  in  the  initial  stage.  This  mainly  occurs 
because  the  initial  amplitude  adn  of  the  0(l)-viscosity  solution  approaches  of  the  viscous-limit 

solution.  The  growth  rates  are  independent  of  A  in  the  initial  linear  stage.  When  A  =  105,  the  finite- 
viscosity  solution  closely  follows  the  viscous-limit  solution  during  the  initial  linear  stage  and,  also,  in 
the  beginning  portion  of  nonlinear  stage,  then  it  jumps  onto  the  asymptotic  solution  rather  abruptly. 
Although  they  are  plotted  for  different  values  of  A,  in  order  to  compare  against  the  corresponding 
finite- viscosity  solutions,  A  is  not  an  independent  parameter  for  the  viscous-limit  solutions. 
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Figure  9.  Inviscid  2  +  2  mode  critical- layer  interaction  when  9P  =  30°,  6S  =  50°  and  A  =  0.  (a) 
In  \Ap\  (black),  ln|As|  (green);  (b)  ln|v4di|  (blue),  In  |Ad2|_(red);  (c)  Re{A!dJ Adl)  (blue),  R e(A'd2/Ad2) 
(red);  (d)  hn.(A'dlf Ad{)  (blue),  lm.(A'd2/ A^)  (red)  versus  t.  Solid  curves,  0(l)-viscosity  solutions;  dot- 
dashed  curves,  initial  exponentially-growing  solutions;  dotted  curves,  asymptotic  solutions.  Blue  and 
red  dotted  curves  are  indistinguishable  and  only  red  dotted  curves  are  visible  in  (c,  d). 


Figure  10.  Effect  of  viscosity  on  a  2  +  2  mode  interaction  when  9P  =  30°  and  6S  =  —57°.  A  =  0  (inviscid, 
black),  10  (red),  103  (blue)  and  105  (green),  (a)  In  | |  and  (b)  In  | ^4^1  versus  t.  Solid  curves,  0(1)- 
viscosity  solutions;  dashed  curves,  viscous-limit  solutions;  dotted  curves,  asymptotic  solutions. 
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Figure  11.  Effect  of  viscosity  on  a  2  +  2  mode  interaction  when  9p  =  30°  and  6S  =  50°.  A  =  0  (inviscid, 
black),  10  (red),  103  (blue)  and  105  (green),  (a)  In  |^4rf! |  and  (b)  In  | 1  versus  t.  Solid  curves, 
0(l)-viscosity  solutions;  dashed  curves,  viscous-limit  solutions;  dotted  curves,  asymptotic  solutions. 


Figure  11  shows  that  large  viscosity  also  delays  the  onset  of  the  nonlinear  growth  when  6P  =  30° 
and  9S  =  50°.  The  primary  and  secondary  amplitudes  are  the  same  as  those  plotted  in  Figure  9(a). 
Note  that  there  are  spikes  in  the  viscous-limit  solutions  plotted  as  dashed  curves  in  Figure  11.  They 
are  due  to  near  singular  behavior  of  the  viscous-limit  solutions  (119).  The  viscous-limit  solution  can 
actually  become  singular  for  certain  combinations  of  9p  and  9S. 

Numerical  results  of  an  inviscid  2  +  3  mode  interaction  are  plotted  in  Figure  12.  When  9P  =  30° 
and  9S  =  —70°,  the  2  +  3  mode  interaction  generates  three  difference  modes  with  propagation  angles 
of  (' Odi ,  9d2,  9d3)  ~  (—78.5°,  —86.2°,  86.9°).  The  primary,  secondary  and  difference  modes  denoted  by 
P ,  S ,  D i,  D'2  and  D3  in  Figures  2(g,  h)  correspond  to  this  specific  case.  As  plotted  in  Figure  12(a), 
the  first-difference  mode,  in  addition  to  the  primary  and  secondary  waves,  grows  exponentially.  The 
last  two  difference-mode  amplitudes,  Ad2  and  Ad3  plotted  as  solid  curves  in  Figure  12(b),  grow  faster 
than  exponentially. 

The  viscous  effects  on  the  above  2  +  3  mode  interaction  are  shown  in  Figure  13.  The  linear  primary 
and  secondary  wave  amplitudes,  which  are  independent  of  A,  are  given  in  Figure  12(a).  Although 
the  growth  rate  of  the  exponentially-growing  first-difference  mode  is  not  affected  by  the  viscosity  as 
shown  by  (97),  its  initial  amplitude  determined  by  (98)  is  a  function  of  A.  The  complex  ad\  becomes 
(-1.7  -  1.3i)  x  1(T3,  (0.2  -  2.2i)  x  10“3,  (9.4  +  O.li)  x  1CT5  and  (1.1  +  0.5i)  x  10"6,  when  A  =  0,  10, 
103  and  10'5,  respectively.  Effects  of  A  on  the  growth  of  the  2nd  and  3rd  difference  modes  are  shown  in 
Figure  13. 

Figure  14  shows  the  numerical  results  of  the  inviscid  2  +  4  mode  interaction  when  9p  =  30°  and 
9S  =  —76°.  The  2  +  4  mode  interaction  generates  four  difference  modes  with  propagation  angles 
of  (9di,  9d2,  9d3,  9di)  ~  (—80.2°,  —84.2°,  —87.9°,  88.6°).  The  linear  primary  (black)  and  secondary 
(green)  waves  and  the  first  (purple)  and  second  (brown)  difference  modes  are  shown  in  Figure  14(a). 
Their  growth  rates  are  /%,(=  1),  fcs(~  2.87  —  13.71i),  Kdl  =  k*  +  kp  and  Kd2  =  Kdl  +  np.  Also,  the 
initial  growth  rates  of  the  third  (blue)  and  fourth  (red)  difference  modes  become  Kd3  =  nd2  +  kp 
and  ndi  =  Kd3  +  kp  as  shown  in  Figure  14(b).  The  amplitudes  Ad,3  and  Ad4  are  first  reduced  at  the 
beginning  of  the  nonlinear  growth  stage,  but  they  soon  approach  the  asymptotic  solutions  plotted  as 
dotted  curves.  Figure  15  shows  the  evolution  of  Ad3  and  Adi  for  various  values  of  A. 

Effects  of  secondary- wave  propagation  angles  on  2  +  77  mode  interactions  are  investigated  by  vary- 
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Figure  12.  Inviscid  2  +  3  mode  critical-layer  interaction  when  9p  =  30°,  9S  =  —70°  and  A  =  0.  (a)  In  \AV\ 
(black),  In  | .As |  (green),  ln|Adi|  (purple);  (b)  ln|Arf2|  (blue),  ln|A.d3|  (red)  versus  t.  Solid  curves,  0(1)- 
viscosity  solutions;  dot-dashed  curves,  initial  exponentially-growing  solutions;  dotted  curves,  asymp¬ 
totic  solutions. 
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Figure  13.  Effect  of  viscosity  on  a  2  +  3  mode  interaction  when  9P  =  30°  and  9S  =  —70°.  A  =  0  (inviscid, 
black),  10  (red),  103  (blue)  and  105  (green),  (a)  In  | A^2 |  and  (b)  ln|Arf3|  versus  t.  Solid  curves,  0(1)- 
viscosity  solutions;  dashed  curves,  viscous-limit  solutions;  dotted  curves,  asymptotic  solutions. 
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Figure  14.  Inviscid  2+4  mode  critical-layer  interaction  when  9p  =  30°,  9S  =  —76°  and  A  =  0.  (a)  In  \A 
(black),  In  | A; |  (green),  In  |  Ail  (purple),  ln|A2|  (brown);  (b)  In  | Ad3 |  (blue),  ln|A4|  (red)  versus 
Solid  curves,  0(l)-viscosity  solutions;  dot-dashed  curves,  initial  exponentially-growing  solutions;  dotted 
curves,  asymptotic  solutions. 


Figure  15.  Effect  of  viscosity  on  a  2+4  mode  interaction  when  9P  =  30°  and  9S  =  —76°.  A  =  0  (inviscid, 
black),  10  (red),  103  (blue)  and  105  (green),  (a)  ln|As|  and  (b)  In  |  Ail  versus  t.  Solid  curves,  0(1)- 
viscosity  solutions;  dashed  curves,  viscous-limit  solutions;  dotted  curves,  asymptotic  solutions. 
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ing  Qs  while  keeping  9P  unchanged.  The  value  of  0S  is  varied  from  —89°  to  89°  with  increments  of  0.05° 
for  the  inviscid  and  0.1°  for  the  viscous  case.  The  series  solutions  are  numerically  evaluated  for  each 
set  of  (dp,  9S)  combinations.  The  results  from  various  sets  of  (9P,  9S )  combinations  are  then  assembled 
to  obtain  the  amplitudes  as  functions  of  0S  at  given  values  of  t.  Figure  16  shows  the  results  of  this  pa¬ 
rameter  study  when  9P  =  30°.  The  absolute  values  of  the  amplitudes,  \At\,  where  i  =  {p,  s,  dl,  d2,  d3} 
for  the  primary,  secondary,  lst-difference,  2nd-difference  and  3rd-difference  modes,  are  plotted  in  the 
first/leftmost  to  fifth/rightmost  columns  of  sub- figures.  The  effects  of  viscosity  are  also  investigated  by 
varying  A  =  0  (inviscid,  lst/top  row),  10  (2nd  row),  103  (3rd  row)  and  105  (4th/bottom  row).  In  the 
figures,  \AZ\  versus  9S  are  plotted  at  t  =  0,  3,  6,  9  and  12  (only  when  A  =  105).  The  wave  propagation 
angles  of  difference  modes,  9<ii,  0d2  and  Odz,  as  functions  of  9S  for  a  constant  9p  =  30°  can  be  found  from 
Figure  6.  For  the  fixed  value  of  9p  =  30°,  the  amplitude  of  the  linear  primary  wave,  Ap,  is  independent 
of  9S  and  A  as  shown  in  the  first  column  of  sub-figures  in  Figure  16. 

The  2  +  1  mode  interaction  occurs  when  0  <  |@s|  <  30°  as  one  can  find  from  Figure  5.  The 
cyan-color-filled  sections  in  the  color  bar  plotted  on  top  of  each  sub-figures  in  Figure  16  correspond 
to  these  ranges  of  9S.  The  2  +  1  mode  interaction  produces  nonlinear  growth  effects  on  the  secondary 
and  (1st)  difference  mode  amplitudes,  As  and  J4^i(=  A ^),  as  shown  in  the  2nd  and  3rd  columns 
of  sub- figures.  Since  the  2  +  1  mode  interaction  only  generates  a  single  (1st)  difference  mode,  the 
corresponding  regions  to  the  cyan-color-filled  sections  in  the  4th  and  5th  columns  of  sub-figures,  i.e. 
for  the  2nd  and  3rd  difference  modes,  are  empty.  Differences  between  the  results  of  the  inviscid  2  +  1 
mode  interactions  in  Figures  16  (b,  c)  and  those  plotted  in  Figures  10(a,  b)  in  Lee  (2012)  are  due  to 
different  normalizations  of  variables.  The  nonlinear  effect  is  stronger  when  9S  is  negative  than  when  it 
is  positive.  The  effect  of  viscosity  is  shown  to  delay  the  onset  of  the  nonlinear  growth.  The  value  of  9S 
for  the  largest  |A,|  and  \Adi\  at  t  =  9  becomes  closer  to  —  9P  as  the  viscous  effect  is  increased.  Details 
of  amplitude  evolution  when  ( 9P ,  6S )  =  (30°,  —15°)  can  be  found  in  Figures  6  and  8  in  Lee  (2012). 

If  30°  <  |#s|  <  64.34°,  that  corresponds  to  the  yellow-color-filled  sections  in  the  color  bars,  the  2  +  2 
mode  interactions  occur.  The  secondary  wave  remains  linear  as  shown  in  the  2nd  column  of  sub-figures. 
Its  linear  growth  rate  that  is  independent  of  A,  can  be  found  in  Figure  1  in  Lee  (2012).  The  amplitudes 
of  the  first  and  second  difference  modes  are  nonlinearly  affected.  The  amplitudes  |A^i|  and  |  A(j2  \  grow 
faster  for  negative  values  of  9S  compared  to  the  positive  ones  as  in  2  +  1  mode  interactions.  The  value 
of  6S  =  —57°  for  the  detailed  study  presented  in  Figures  8  and  10  was  chosen  to  show  the  results  near 
the  strongest  2  +  2  mode  interaction  when  9P  =  30°. 

The  2  +  3  mode  interaction  occurs  when  64.34°  <  |0S|  <  73.22°,  or  in  the  orange-color-filled  sections 
in  the  color  bars  in  Figure  16.  Both  secondary  and  first-difference  modes  grow  exponentially  as  shown 
in  the  2nd  and  3rd  columns  of  sub- figures.  The  amplitudes  of  the  second  and  third  difference  modes, 
plotted  in  the  4th  and  5th  columns,  are  nonlinearly  affected. 

If  73.22°  <  |0S|  <  77.5°,  as  in  the  black-color-filled  sections,  the  2  +  4  mode  interaction  occurs. 
The  amplitudes  As,  A^i  and  A$2.  grow  exponentially.  The  nonlinear ly-growing  A dz  is  plotted  in  the 
rightmost  column  of  sub-figures.  The  fourth-difference  mode  amplitude,  A is  not  shown  in  the  figure. 

There  are  special  combinations  of  ( 9p ,  0S )  for  which  the  nonlinear  effects  are  reduced  or  disappear. 
When  0S  =  0,  the  linear  plane  secondary  wave,  that  propagates  along  the  downwind  direction,  grows 
exponentially.  The  difference  mode  is  not  generated  by  the  nonlinear  interaction  of  the  oblique  primary 
and  plane  secondary  waves.  If  9S  =  —9P,  the  linear  secondary  mode  grows  exponentially  as  shown  in 
Figure  16.  The  difference  mode  becomes  steady  and  spanwise  periodic.  When  0S  =  9P,  the  nonlinear 
interaction  of  the  primary  oblique  wave  interacting  with  itself  does  not  generate  a  difference  mode 
(Goldstein  &  Leib  1989).  If  9S  equals  the  plus/minus  wave  propagation  angle  corresponding  to  as  = 
Nap  for  N  >  2,  the  first  IV  —  1  difference  modes  grow  exponentially  as  confirmed  in  Figure  16.  The 
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IVth-difference  mode  becomes  steady  and  spanwise  periodic. 

It  is  found  that  a  2  +  2  mode  interaction  can  generate  a  plane  2nd-difference  mode  when  the  relation 
(131)  is  satisfied.  With  the  parameter  values  in  (14),  the  pair  of  ( 9P ,  9S)  ~  (30°,  43.87°)  satisfies  the 
condition  independently  of  A. 

In  Figure  7,  we  find  that  both  coefficients  < pa  nv^  and  (pAnv'>  in  the  inviscid  amplitude  equations  (109) 
and  (110)  vanish  when  the  2nd-difference  mode  is  a  plane  wave,  the  case  already  considered  above. 
Figure  7  also  shows  that  <$>£nv'>  =  0  but  (p^nL^  /  0,  for  example,  if  ( 9P ,  9S )  equals  (30°,  54.21°)  for  a 
2  +  2  mode  interaction.  The  amplitudes  A^i  and  grow  exponentially  as  shown  in  Figures  16(c,  d). 
We  also  find  from  Figure  16  that  this  inviscid  effect  disappears  as  A  is  increased. 

In  Figures  17  and  18,  the  results  of  a  parameter  study  when  9p  =  50°  and  70°,  respectively,  are 
summarized.  Effects  of  9S  on  various  2  + IV  mode  interactions  are  investigated  by  varying  9S  from  —89° 
to  89°  as  in  the  previous  9P  =  30°  case.  The  regions  of  various  2  +  N  mode  interactions  are  shown  as 
color-filled  sections  in  the  color  bars.  The  values  of  1 6S  |  at  the  boundaries  where  as  =  Nap  become 
50°,  71.25°,  77.63°  and  80.75°  if  6P  =  50°  for  IV  =  1,  2,  3  and  4  as  can  be  found  from  Figure  5.  They 
are  70°,  80.15°,  83.45°  and  85.09°  if  9p  =  70°.  Differences  between  the  results  of  the  inviscid  2  +  1 
mode  interactions  in  Figures  18  (b,  c)  and  those  plotted  in  Figures  11  (a,  b)  in  Lee  (2012)  are  due  to 
different  normalizations  of  variables.  The  results  of  these  9P  =  50°  and  70°  cases  are  consistent  with 
those  of  the  9P  =  30°  case.  The  nonlinear  parametric  growth  is  more  effectively  enhanced  when  the 
sign  of  the  propagation  angle  of  the  secondary  wind  wave  9S  is  opposite  to  that  of  the  primary  wind 
wave  9p  than  when  both  have  the  same  sign. 


38 


\n\Ap\vs.6a  ln|As|  vs.  9S  ln|Adi |  vs.  9S  ln|Ad2|  vs.  6S  ln|A<H|  vs. 
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\n\Ap\vs.6a  ln|As|  vs.  9S  ln|v4di  |  vs.  6S  ln| Ad2 1  vs.  6S  lnlA^I  vs. 
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\n\Ap\vs.6a  ln|As|  vs.  9a  ln|j4di  |  vs.  9S  ln|A,a|  vs.  9a  ln|A<H|  vs. 
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DISCUSSION  AND  CONCLUSIONS 


A  nonlinear  interaction  that  can  generate  low-frequency  fluctuations  in  the  air  above  obliquely 
propagating  wind  waves  is  studied.  The  nonlinear  interaction  can  be  initiated  by  two  three-dimensional 
fluctuations  in  wind,  which  are  induced  by  wind-driven  surface  waves,  referred  to  as  the  primary  and 
secondary  waves.  The  secondary-wave  frequency  can  be  higher  than  the  primary  frequency.  The 
analysis  is  valid  for  general  two-dimensional  wind  profiles.  The  characteristic  wave  speed  is  assumed  to 
be  sufficiently  large,  so  that  the  critical  layer  is  distinct  from  the  viscous  sublayer.  The  experimental 
study  by  Hristov,  Miller  &  Friehe  (2003)  has  shown  the  existence  of  distinct  critical  layers  in  an  open 
ocean. 

As  pointed  out  in  Lee  (2012),  the  phase  speed  in  the  wave  propagation  direction  of  an  oblique  wave, 
A,  becomes 

Ci  =  uJi/li  =  c  cos  9U  (136) 

for  i  =  {p,  s,  dl, . . . ,  dN},  where  'yl,  6t  and  uj1  are  defined  by  (12),  (13)  and  (32)-(34),  and  c  is  the  wave 
speed  in  the  streamwise  or  downwind  direction  (Drazin  &  Reid  2004,  p.  128).  For  wind  waves  and 
nonlinear  difference  modes  with  an  equal  wave  speed  in  downwind  direction  c,  the  phase  speeds  ct  are 
proportional  to  cos^.  For  an  equal  phase  speed  ct,  the  wave  speed  in  downwind  direction  c  becomes 
larger  and  the  critical  level  moves  higher  from  the  interface  as  the  wave  propagation  angle  is  increased. 
An  oblique  wind  wave  may  have  a  distinct  critical  layer  from  the  viscous  sublayer  although  its  phase 
speed  Ci  is  small. 

We  suppose  that  the  primary  wind  wave  is  the  largest  oblique  wave  in  an  initial  wind-driven  wave 
field.  The  linear  primary  wave  grows  exponentially.  According  to  Miles’  (1957)  theory,  its  growth 
rate  is  proportional  to  a.  that  is  the  density  ratio  of  air  to  water  as  in  (1).  The  typical  value  of  a 
is  0(1CR3).  Once  the  steepness  of  the  primary  wind  wave  becomes  0(<r3),  it  can  nonlinearly  interact 
with  a  secondary  oblique  wave  of  smaller  steepness.  If  the  wave  speeds  in  downwind  direction  of  the 
primary  and  secondary  oblique  waves  are  equal,  both  wind  waves  share  a  common  critical  layer  where 
the  nonlinear  interaction  first  takes  place  (Davis  1969;  Benney  &;  Bergeron  1969;  Drazin  &;  Reid  2004). 
If  N  is  the  smallest  integer  not  less  than  the  ratio  of  the  secondary  to  primary  frequencies,  or  if  the 
frequency  ratio  is  greater  than  N  —  1  but  less  than  N,  as  in  (15)  and  (16),  the  critical-layer  interaction 
generates  N  difference  modes.  The  2  +  N  mode  interaction  is  composed  of  two  wind  waves  and  N 
nonlinearly  generated  difference  modes.  A  plane  surface  wave,  propagating  in  the  streamwise  direction, 
whose  steepness  is  larger  than  the  primary  oblique  wave  can  coexist  without  playing  any  significant 
role  in  the  present  nonlinear  interaction  stage. 

When  N  >  2,  the  interaction  between  the  primary  and  secondary  wind  waves  generates  the 
first-difference  mode.  The  primary  wave  then  couples  with  the  first- difference  mode  to  generate  the 
second-difference  mode.  This  sequential  difference-mode-generation  process  is  continued  until  the  JVth- 
difference  mode  is  generated  by  the  quadratic  interaction  of  the  primary  and  (AT—  l)th-difference  modes. 
The  critical-layer  interaction  between  the  primary  and  iVth-difference  modes  produces  a  nonlinear  ef¬ 
fect  on  the  (N  —  l)th-difference  mode.  The  amplitudes  of  difference  modes  are  of  the  same  order  of 
magnitude  as  the  secondary-wave  amplitude  in  air.  Both  primary  and  secondary  wind  waves  remain 
linear.  The  first  N  —  2  difference  modes,  whose  frequencies  are  higher  than  the  primary  frequency  as  in 
(36),  grow  exponentially  although  they  are  nonlinearly  generated.  The  growth  rate  of  the  nth-difference 
mode,  or  initial  growth  rate  of  the  (N  —  l)th  or  Nth  difference  mode,  is  increased  by  nnp  from  k*  as 
shown  in  (97)  for  1  <  n  <  N .  Meanwhile,  its  frequency  u(jn  is  decreased  by  nojp  from  the  highest  us 
as  obtained  from  (34)  for  1  <  n  <  N  —  1.  The  np  and  ks  are  the  linear  growth  rates  of  the  primary 
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and  secondary  waves,  and  ojp  (=  avc )  and  ujs  (=  asc)  are  their  frequencies.  The  amplitudes  of  the 
(TV  —  l)th  and  TVth  difference  modes  are  governed  by  coupled  integral  equations,  if  they  do  not  satisfy 
the  free-surface-wave  dispersion  relation  (11).  Numerical  results  show  that  they  can  grow  faster  than 
exponentially.  In  fact,  both  of  them  can  grow  like  the  exponential-of-an-exponential  in  the  inviscid 
and  finite- viscosity  cases.  The  frequencies  of  the  last  two  difference  modes  are  lower  than  the  primary 
frequency  as  shown  by  (36).  In  the  special  case  where  the  relation  (11)  is  satisfied  by  a  difference  mode, 
the  amplitude  will  be  governed  by  an  integro-differential  equation  that  is  similar  to  that  of  a  secondary 
wave  in  Lee  (2012). 

The  nonlinear  growth  in  a  2  +  TV  mode  interaction  is  of  the  parametric-growth  type.  The  coupling 
between  the  primary  and  TVth-difference  modes  produces  a  nonlinear  effect  on  the  (TV  —  l)th-difference 
mode  if  TV  >  2.  When  TV  =  1,  the  quadratic  interaction  between  the  primary  and  difference  modes 
produces  a  nonlinear  growth  of  the  secondary  mode. 

It  is  shown  that  the  2  +  TV  mode  critical-layer  interaction  greatly  enhances  the  growth  of  oblique 
modes  whose  frequencies  are  lower  than  the  primary  frequency.  They  are  the  (TV  —  l)th  and  TVth 
difference  modes  if  TV  >  2,  and  the  secondary  wind  wave  and  difference  mode  when  TV  =  1.  Due 
to  nonlinear  effects,  the  amplitudes  of  these  modes  can  quickly  become  as  large  as  the  primary-wave 
amplitude,  i.e.  0(a'i).  Their  nonlinear  coupling  in  a  following  stage  can  generate  multiple  difference 
modes  whose  amplitudes  are  0(ai).  The  nonlinear  interaction  in  this  stage  will  still  be  of  a  quadratic- 
interaction  type,  by  which  any  combination  of  two  modes  can  interact  to  generate  a  difference  mode 
or  produce  a  nonlinear  effect  on  it.  The  most  amplified  modes  whose  amplitudes  first  become  0(a 5/2) 
may  then  enter  the  next  self-interaction  stage.  Lee  &  Wundrow  (2011)  show  that  these  modes  in  the 
self-interaction  stage  can  grow  even  faster  than  they  do  in  the  parametric-growth  stages. 

Although  the  2  +  TV  mode  interaction  is  induced  by  wind-driven  surface  waves,  the  most  rapidly 
growing  mode  in  air  may  not  be  the  one  that  is  associated  with  a  wind  wave.  For  example,  if  TV  >  2, 
the  (TV  —  l)th  and  TVth  difference  modes,  that  can  grow  faster  than  exponentially,  do  not  satisfy  the 
free-surface-wave  dispersion  relation.  The  amplitudes  of  the  corresponding  disturbances  in  water  are 
reduced  by  a  factor  of  0(a)  as  shown  in  (28)  and  (42).  For  this  interaction,  there  might  be  hardly  any 
observable  evidence  of  nonlinear  effects  on  the  water  surface  although  the  difference  modes  in  air  have 
already  become  as  large  as  the  primary  wave.  This  indicates  that  the  multi-mode-generation  process 
in  air  may  start  earlier  than  in  water. 

The  critical-layer  interaction  could  play  an  important  role  in  the  growth  of  perturbations  in  the 
wind  over  water.  One  of  the  results  of  the  present  analysis  is  that  the  2  +  TV  mode  critical-layer 
interaction,  occurring  in  the  air  over  wind-driven  surface  waves,  can  nonlinearly  enhance  the  growth  of 
air  fluctuations  whose  frequencies  are  lower  than  the  primary  frequency. 
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APPENDIX  A:  RELATIONS  BETWEEN  Jz  AND  bt 

For  i  =  {p,s,dl, . . .  ,dN},  <3?*  can  be  computed  by  numerically  integrating  (44)  subject  to  (46)  and 
(49)  below  the  critical  level  where  0  <  z  <  zc.  Then  the  complex  constant  Jt  defined  by  (48)  can  also 
be  evaluated  as 

Jt  =  (1  /6,)  lim  \z ).  (A-l) 

Z->ZC 

The  second  shape  function  becomes  purely  imaginary  on  the  real  axis  below  the  critical  level 
since  the  boundary  conditions  (46)  and  (49)  are  purely  imaginary.  Thus,  JA  also  becomes  purely 
imaginary,  or 

Re(  JA)  =  0.  (A-2) 

By  multiplying  (44)  by  integrating  over  transverse  domain  with  (52),  and  substituting  (86),  we 

obtain 

£  \dMb)\2  +  £  +  ££j  i^(6)i2  dz  =  -j:k  £ja£  -  £ .  (a-3) 

Because  the  left-hand  side  of  the  above  equation  is  real,  the  term  in  the  parentheses  on  the  right-hand 
side  must  be  purely  imaginary.  Therefore, 

Ke^nJAU” /U'c  -  1/6,)  =  0  and  \Jt\  =  |t//Re(6,)/(7r[/"6,3)|.  (A-4) 

If 

U"/(U  -  c)  +  7t2  >  0  (A-5) 

below  the  critical  level,  which  is  satisfied  by  the  logarithmic  wind  profile  (7),  then 

d2Mb)/^b)  >  0.  (A-6) 

From  above  results  with  (A-l)  and  (A-6),  one  can  show  that  (A-2)  becomes 


Jib  i  —  ]\Jib% 


(A-7) 


By  putting  (A-7)  into  (A-4),  we  find  that 

Jt  =  -it//Re(6,)/(7rt/"|6i|26,). 


(A-8) 
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APPENDIX  B:  SYSTEM  OF  CRITICAL-LAYER  EQUATIONS 


The  normalized  variables  in  (79)-(81)  and 

{a,  qj  =  {A/*,  7?4(1) 


(B-l) 


(2)  . 


are  used  for  the  critical-layer  equations.  The  variables  represented  by  Q3  '  in  (78)  are  also  normalized 
as 

{cff>,rw,iyf,pf}  =  {2t/f,  vf ,  ±w?\  (b-2) 

with  dj  given  below  by  (B-7). 

The  system  of  normalized  equations,  which  we  refer  to  as  the  critical-layer  equations,  now  becomes 
at  the  leading  order,  for  i  =  {p,  s,  dl, . . . ,  dN}, 


—  d^A%, 

subject  to  the  transverse  boundary  condition,  Ql~l/(>  as  »oo.  We  have  put 

„  d  ...  x  d2 

Cl  =  -+iai(-\^. 

The  second-order  critical-layer  equations  are  written  as 

£j{V}2\  W%1  uW}  =  {R®,  R$,  R$j, 

WH  =  &3U32)  +  fa™  >  PH  =  °> 
for  j  =  {p,  0,  2 p,  s,  p  +  s,  dl,  . . . ,  d(N+ 1) } ,  with 

bo  —  0,  «2P  —  2bp,  bp_(_s  —  dtp  T  ds,  otdn  —  Tidp  bs, 

Po  =  0,  /?2; p  =  2/3p,  Pp+s  =  Pp  T  Ps  f  Pdn  =  Tlpp  —  Ps. 

For  the  primary,  secondary  and  difference  modes  (B-5)  and  (B-6)  become 

CpV^  =  -i  pppj?)  -  \{a/P)^MjWc  U'P  dpPp{  C  +  CofQp , 

£pWW  =  i(b/ P)y/ M/U'c  UPdpPpAp, 

w$  =  dpuM  +  ppvW ,  PjV  =  0, 

CsvP  =  -i PsPP  -  \{d/P)^/M/U'c  U" as£is((  +  lofQs 

—ppPdlApQdlc  ~~  PdlPp-^dlQpc  ~  d^d2ppPd\QpQd\i 


(B-3) 

(B-4) 

(B-5) 

(B-6) 

(B-7) 

(B-8) 

(B-9) 

(B-10) 

(B-ll) 

(B-12) 


L  cc 


CsWjV  =  i(d/ P)^MjU>  up  as% As  -  ^  [dp%pdlApQ*dl 

Add\pdlPp-^dlQp qq  T  ^^p2d^p^dlPpPdl{QpQdl)c 


(B-13) 


B-l 


wW=asUP  +  psVsM,  Pjf=  0 


(B-14) 


£dnVdn  =  -i PdnPdn  ~  \  («/A)  VM/Uc  Uc  «dnAdn(C  +  Co)2<3dn  -  >Ad(n+l)^p<2d(n+l)  ( 

7d(»i+l)Ap74d(r)+l)Qp^  iL'p—d  &d(n-\-2)V,pH,d(n-\-l)QpQd(n+l)  ~^‘PpP'd(n—l)-^-pQd(n—l)( 
^~'Jd(n—l)P'p-^d(n—l)Qpt;  i  ^p—d  ^d(n—2)Ppfdd{n—l)QpQd{n—l)i 


(B-15) 


^dnWdn^  =  i(®/A)  VM/U'  Uc  Qtdn'JdnAdn  &p_d  ^  &plfpP'd(n+l)ApQd(n+l)  ^ 
_^®d(n+l)7d(?i+l)Ap^d(n+l)Qp(-(-  l^V p_d  ^p^d(n+l)^'P^d(n+l)^QpQd(n+l)')C 

~^-^p—d  *-^pHp^d(n— l)-ApQd(n— 1)  ^  ®c£(n—  1)7(2(71— l)Ap^-d(n—  l)Qp(( 

~l^p_d  ^ ®P®d(n—  l)ApAd(n— l)(QpQd(n—  l))c  > 

w(2)  _  ^  7j(2)  i  a  t/(2)  p(2)  _  n 

VVdnc  ~  adn^dn  +  PdnV dn  ,  rdnc  ~  U’ 

where 

A*  =  A/7*  =  sin  dt  for  *  =  {p,  s},  pdn  =  fidn/ldn  for  n  =  1,  2, . . . ,  N  +  1, 

I'p—d  =  ftp/ Pdn/^dn- 


(B-16) 

(B-17) 

(B-18) 

(B-19) 
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APPENDIX  C:  NORMALIZATION  PARAMETERS  a,  (3,  k,  M  AND  Co 

The  normalization  parameters  are  chosen  as 

«  =  «30>  P  =  6l30  =  7r/6,  (0-1) 

k  =  — 7rcC/C|^3o|2(cos 930)/[2aU'2  tanh(730-ff)],  (C-2) 

Co  =  -U^R3o/(irU,’a30\bao\2),  VM  =  P2\ J30||630|2(cos@30)/(a2v/^),  (C-3) 

where  the  subscript  30  denotes  the  quantities  of  a  surface  wave  with  propagation  angle  of  30  degrees. 
The  choice  of  #30  =  7r/6  is  arbitrary.  According  to  Fjprtoft’s  theorem,  it  is  assumed  that  U"  < 
0  and  U'c  >  0,  so  that  k  and  \[M  are  positive.  From  (79)  and  (C-l),  one  can  show  that 

®30  =  ^30  =  1-  (C-4) 
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APPENDIX  D:  CRITICAL-LAYER  JUMPS  FOR  THE  DIFFERENCE 

MODES 


If  the  solution  (89)  is  substituted  for  Ap,  the  variables  A^(t),  F±(t)  and  G^(t)  in  (83)  become 


An  =  ±  iir(a/ (3) \J M/U'c  U"(c(dn/6idn)Adn(t), 


(D-l) 


Fn  =  =F  Gp  Unld{n±  1)  /o°°  Ea±  (t)  (j)  T2  Ad{n±1)  (t  -  t)  dr 

-Cfe  fo°°  r  En±  (n,  T2)  e£L\t!  +  r2)  n  Ad(n±1)(f  -  n  -  r2)  * 


(D-2) 


G«  =  =F 


W 


where 


Vnlp  r  Eb±(T)  Ec±\r)  T2  Ad(n±1)(f  -  r)  dr 


r  fo°  Ev±  iTuT2)  Ec'1\ti  +  r2)  n  Ad(n±1)(f  -  n  -  r2)  dndr2 


un  =  n-v ,  i>  =  as/ap 


Ca±  =  2tt(1  ±  Vn)3a5l'p!51)Ad(n±l).  Cfe  =  2ad(n±i)z/^1)/ip, 

=  7rapc£},  C™  =  ±  z/n)2], 


(D-3) 

(D-4) 

(D-5) 


*(r)  =  exp  [-|Aa2r'2(l  ±  z'n/r3 


’(r )  =  exp  [~\\a2v2{l  ±  un)2TZ 


^(r)  = 


(r)  =  exp  [-(1  ±  t'n)^r 


£|i±) (ri,  r2)  =  exp  [-gAd2  {z^rf  +  (1  ±  z^rf  +  {±z/„n  +  (1  ±  z^n)r2}3}]  , 
Ei±{n,T2)  =  exp  [— |Aa2(l  ±  z/„)2  {(1  ±  vn)(VnTi  +  Tl)  +  (T^nTi  +  r2)3}]  , 
with  dp  =ap,  d^  =  a*,  d+  =  /%,,  =k ;*,  Arf0  =  -A*  and  /ip,  //<*„,  vpr^d  in  (B-18)  and  (B-19). 


(D-6) 

(D-7) 

(D-8) 
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APPENDIX  E:  VARIABLES  gn-i,  hn,  sasy  AND  rasy 

When  A  is  0(1),  gn-\  and  hn  in  (124)  and  (125)  become 


9n- 1  — 


16(R e/ip)6(n  -  1  +  £)3(n  -  1  +  £  +  vN)3  [h{n)  +  h{n)\  [i|(n-l)  +  7|(n-l)] 


^(jvWjv-1) 

'-'a- 


"*7?  -  C(bN) 


UN- i7p  -  C'h  0 


=  4(Re«p)3(n  +  0"  [4(n)  +  I4(n)]  /  C{aN_\vN%  -  C£’ / 2) 


'•(n) 


where 


A(n)  =  r  ^(n)  f2dfu 

I2{n)  =  <£*<£*  /0°°  f0°°  4+_l)(n,  r2)  e-[2(n+iV-1)Re«p-^;+«3*](f1+f2)  fi  d^dflj 
/3(n)  =  C[NJun%  f1)e-(2nReKp+z>«p+«s*)f1  f2dflj 

I4(n)  =  -C^Jc^  f0°°  /0°°  4-}(ri,  f2)  e-(2nReKp+^+^‘)(fi+f2)  f4  df2dfi, 
and  £,  0,  un,  C%±  ,  Cb'^\  and  are  given  by  (126),  (D-4)-(D-6)  and  (D-8). 


The  coefficients  sasy  and  rasv  in  (129)  and  (130)  are 


asy 


Sa,sy 


S0/(8V3  7 r5/2)]  r3(C)r3(e  +  uN)  {(^/Re«p)3^)/4} 


-£,-vN/2+7/12 


3  fTb — 1  ^ 


m=l 


V 

asy 


vt") 


~(inv) 

ra 


-1/2 


(E-l) 

(E-2) 

(E-3) 

(E-4) 

(E-5) 

(E-6) 


(E-7) 

(E-8) 


E-l 
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